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Abstract 

Anomalies and BRST invariance are governed, in the context of Lagrangian Batalin-Vilkovisky 
quantization, by the master equation, whose classical limit is (S, S) — 0. Using Zimmerman's 
normal products and the BPHZ renormalisation method, we obtain a corresponding local 
quantum operator equation, which is valid to all orders in perturbation theory. The formu- 
lation implies a calculational method for anomalies that is useful also outside this context 
and that remains completely within regularised perturbation theory. It makes no difference 
in principle whether the anomaly appears at one loop or at higher loops. The method is 
illustrated by computing the one- and two-loop anomalies in chiral W3 gravity. 
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1 Introduction 



The most comprehensive method for treating gauge theories is at present the method of 
Batalin and Vilkovisky (BV). It unifies in one framework the description of gauge invariance, 
the selection of physical observables, the gauge fixing, BRS methods, the introduction of ghosts 
(and antighosts, if wanted), etc. All this is achieved irrespective of whether the gauge algebra 
is a relatively simple one (as for nonabelian gauge theories) or a more complicated one, closing 
only upon using the equations of motion (a BRST charge that is nilpotent only on shell), as 
often occurs in supergravity and in string theories. 

Whereas over the years the quantum BRS treatment of gauge theories |Qj has met with many 
successes !|, the quantum theory in the BV scheme is less developed. Indeed, the remarks 
above strictly apply to the classical theory. For many considerations, including some quantum 
questions, like the classifications of one loop anomalies for example, one practically needs use only 
of the classical theory: interesting results follow by considering that in the quantum theory the 
anomalies have to satisfy a certain condition and by solving its classical version (the Wess-Zumino 
consistency condition). A full quantum analysis however requires a more careful treatment. It is 
rather easy to give formal arguments, and in fact they were already given in the founding papers 
on the subject Q; but when one goes deeper into the matter, taking seriously the occurrence 
of the well-known divergences, one realises that the literature is defective. This paper addresses 
this problem head-on. 

In the BV setup, the gauge invariance at the classical level is expressed by the so-called 
master equation, (S, S) = 0. Here S is the extended classical action, which consists of the 
classical action (including ghosts and antighosts), source terms for the BRS-transformation 
laws, and some supplementary terms if the gauge algebra is open. On the quantum level, 
BRST invariance can be transcribed as (r, T) = 0, where V is the generating functional of 
the one-particle irreducible graphs, or the effective action, evaluated in the presence of the 
above mentioned source terms. Anomalies consist then in a violation of this quantum equation. 
In the BRS context, this question was developed making extensive use of the formulation of 
the quantum theory on the level of the effective action. This development can of course be 
transcribed to the BV formalism. However, the question arises also how this translates in the 
quantum theory to the local level, i.e. to the level of the action. Formal manipulations modify the 
classical master equation into (S q ,S q ) = 2ifiAS q , where the quantum action possibly contains 
terms that are explicitly of higher order in h, S q = S + 0(K) and A is a second order functional 
differential operator where the functional derivatives are at the same space-time point. If S q is a 
local action, this operation leads to 5(0), it is ill-defined. The source of the problem is of course 
the same as always in local quantum field theory, namely the coincidence of space-time points. 
A regularisation scheme should solve it. 

Efforts in this direction include @ where the method of Pauli and Villars has been applied 
to this question. This proposal works well on the one loop level, but it is not very clear how to 
extend this type of calculation to higher loops. In a less conventional approach, the use of non- 
local regularisation was suggested || g|, providing a method valid also for higher order loops. 
In Q, an analysis of the quantum aspects of BV using dimensional regularization is presented. 

The interest in a full local quantum treatment is not a purely conceptual one. It is also of 
importance for the question of anomalies appearing in higher loops. Indeed, in some theories, 
contrary to the situation in nonabelian gauge theories ||, anomalies can arise at the multiple 
loop level. This is notably the case in extended theories of two dimensional gravity, so-called 
VF-gravities. The presence of a two-loop anomaly in W3-gravity, for example, is ascertained 
through the use of the operator product expansions of two dimensional conformal field theory. 
Whereas this is certainly a valid method, it leaves a gap. Namely, rather than deriving these 
operator product expansions from a proper treatment of the renormalised quantum field theory 
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that underlies them, their validity is mostly assumed, based on algebraic and symmetry consid- 
erations. Clearly, a complete treatment of higher loop anomalies in renormalised perturbative 
quantum field theory is lacking. 

The purpose of the present paper is to formulate and derive the basic local equations govern- 
ing the BV method on the level of renormalised perturbation theory. On a more practical level, 
we will provide a proper method, always within conventional renormalised perturbation theory, 
to compute higher loop anomalies. We will use the method of Bogoliubov, Parasiuk, Hepp and 
Zimmerman (BPHZ) to achieve this purpose. Not only is this method of impeccable standing, 
but also we will find that it is in fact well-suited to put the quantum master equation of the BV 
formalism on a solid foundation. In fact, that it is suited for computing anomalies (in global 
symmetries) was already realised long ago (see ||). Treatments of non-abelian gauge theories, 
using the same operator methods, on the level of the effective action, have of course played an 
important role in firmly establishing their renormalisability (after the dimensional regularisation 
treatment of flO|| ). An overview of this and other developments that goes together well with the 
BV scheme can be found in 

This paper should not be read as a complete account of the way the BV formalism and the 
BPHZ method work together to provide a full account, in renormalised quantum field theory, 
of all possible higher loop anomalies. Rather it should be viewed as the onset of a program. It 
does provide the basic equations and explicit examples. Working out the detailed consequences 
in the most general situations is left for the future. As far as the method to compute BRST 
anomalies is concerned, it can be stripped of its BV dressing and applied outside this context 
as well. 

We have organized the paper as follows. In section || we recall the basics of the Batalin- 
Vilkovisky method, the relationship of the quantum master equation with the occurrence of 
anomalies, and the Zinn-Justin equation, thus setting the stage for a more precise formulation 
of what we want to achieve. Section |3| starts with an extensive discussion of renormalisation 
schemes in the BV setting, thus argueing for the use of the BPHZ method. The basics of this 
method, including the manifestation of anomalies (in particular the chiral anomaly) are then 
recalled, and some aspects of the massless case are presented from an angle that complements 
the existing literature. Section || combines the previous developments into a general derivation 
of the master equation in the BPHZ framework. The corresponding treatment for the concrete 
case of V^Vgravity in section || serves as an illustration of the more abstract treatment of the 
previous sections, and also as a warm-up for the treatment of W3 gravity in section || wherein, 
after quickly giving the one-loop anomaly, we come to the complete computation of the two- 
loop anomaly of W3. Three appendices collect some formulas on BPHZ renormalized one- loop 
integrals for 2D chiral gravities, supply the complete one-loop anomaly for W3, and illustrate 
the BPHZ treatment of the ghost number anomaly. 



2 Batalin-Vilkovisky Lagrangian quantisation 



In this section, we recall the vocabulary of the Batalin-Vilkovisky (BV) formulation of gauge 
theories ||, including the formulation of anomalies on a formal quantum level. More detailed 
reviews on the BV scheme, reflecting the point of view of the present authors, can be found in 
[12, 13, 14]. Other recommended reading on BV is [15]. 



1 This book is also useful for the reader who needs more background material on renormalisation and symmetries 
in field theory than can be provided in this paper. 
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2.1 Basic concepts 



The BV scheme combines BRS quantisation |T6| and its extension to open algebras [17] with 
the observation -originally due to Zinn- Justin (see e.g. chapter 19 of |lq])- that it is convenient 
to introduce sources for the BRST transformation rules during the renormalization process. 
With every fieldQ & A one associates an antifield & A of opposite Grassmann parity^]. Fields and 
antifields are canonically conjugated w.r.t. the antibracket defined by 

1 ' ' 5<S> A 5$ A 5<S> A 5$ A ' { ' 

The ghost number of a field and that of its antifield add up to —1. 

The BRST invariant, classical theory is then described by an extended action 5($,$*), 
formed by adding several pieces: the classical action, a term linear in the antifields Q* A with 
coefficients R A equal to the BRS variation of the corresponding field 

S[$, $*] = 5 [$] + ®* A R A [$\ + . . . , (2.2) 

and additional terms, at least quadratic in the antifields, which encode special features of the 
BRST algebra, e.g. its off-shell structure. This action is to be constructed in such a way that 
the classical master equation 

(S,S) = 0, (2.3) 

is satisfied. A solution of this equation (together with some suitable boundary conditions) always 
exists |Jl9|1 , whereas its local character, universally taken as one of the cornerstones of quantum 
field theory and therefore a desirable characteristic, can always be ensured under reasonable 
conditions on the gauge algebra generators [^0| (see also |2lJ and references therein). 

The antibracket, eq.(2.1), the extended action, eq.(2.2), and the classical master equation, 



eq.([2,3|), play a central role in the theory. For example, the antibracket is used to perform 
canonical transformations of field variables and to implement gauge fixing (which becomes merely 
a change of basis). The extended action is used to define the BRS operation, and the master 
equation (^) guarantees then its on-shell nilpotence. 

2.2 The quantum master equation and anomalies 



In the quantum theory, equation ( |2.3| ) is modified. Indeed, from a purely formal point of 
view, gauge (or BRS) invariance at quantum level, with a given action S, is transcribed into the 
BV formalism as 

±(S,S)=ihAS, (2.4) 



with 



A = mw A ■ (2 - 5 > 



We call attention once again to the implicit notation (see the footnote on page ||) which implies 
that the sum over A in ( |2.5| ) includes not only a summation over all (components of) the fields, 
but also an integration over the space-time point. This point is the same for both functional 



2 $ A is the complete set of fields in the theory: gauge fields, matter fields, ghosts, antighosts, . . . 

3 De Witt notation is assumed throughout the paper: the index A is used to indicate the different fields, their 
components, and the space-time point on which they depend (unless it is explicitly displayed). In this way, a 
summation over A includes not only discrete summations, but also integration over space-time. The derivatives 
are left and right functional derivatives 
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derivatives. Accordingly, the expression eq.( |2.5D is meaningless when acting on a local functional 
-in particular its use in eq. (|2.4j )! In fact, it could be viewed as an expression involving <5(0), 
and one may be tempted to neglect it for that reason, for example appealing to dimensional 
regularisation @|. However, careful consideration leads one to conclude that in fact this term 
is potentially a source of anomalies [fl. In the Pauli-Villars scheme for example, to order h, 
a definite expression can be associated to it. Then one has to decide whether this anomaly is 
genuine: if it can be avoided by adding an extra counterterm %M\, the classical action S should 
be modified to a 'quantum action' S q = S + %M\ + . . . Without restrictions on M\ this is in fact 
always possible (see, for instance, ||), but here the locality requirement comes into play again: 
for local quantum field theory one will refuse to add anything but a local M\ . If one can not find 



a local Mi, the theory has a genuine anomaly, and (2.4) can not be satisfied. This reasoning 



can be extended to higher orders in h. The formal quantum master equation is then 

±(S q ,S q ) = ihAS q , (2.6) 

where S q is the previously mentioned quantum action^, formally reducing to the classical action 
for h — * 0. 

In conclusion, within this framework it is the task of one's regularisation scheme to provide 
a well defined replacement for the purely formal equation (p.6j). 

2.3 The Zinn-Justin equation 

In the quantum theory, the master equation can also be discussed at the level of the ef- 
fective action T[tp, <&*], a functional of the classical fields tp A and of the sources for the BRST 
transformations This functional is obtained from the generating functional of the connected 
diagrams by trading the sources J a for the classical fields ip A by a Legendre transformation, 
while keeping the BRST sources & A fixed. Since T incorporates all quantum corrections to 
the classical action into effective interactions, i.e. T[tp, $*] = S(ip, <£*) + hTi[ip,^*] + the 
quantum BRST structure is naturally described by the equation 

of which the classical master equation can be considered to be the Ti = limit. 



Fulfilment of (2/7) ensures in fact gauge independence of the quantum theory. Instead, if the 



BRST symmetry is broken by loop corrections, i.e. the gauge symmetry is anomalous, equation 



(|2.7[) is replaced, according to Lam's quantum action principle |22j, |23[, by J24J 



\{T,T) = -iK{A-T). (2.8) 

On the right hand side, (.A • V) is the generating functional for 1PI diagrams with one insertion 
of a local operator A[ip, $*]. The BV formalism gives an explicit albeit formal expression for it, 
related to the quantum master equation (fO|), namely A = AS + (i/2Th)(S, S). In this way the 
operator A[(f, <£*] makes explicit the breaking of the classical BRST structure due to quantum 
corrections and is interpreted as the BRST anomaly. The equation (^]7]), in a slightly different 
guise without the sources corresponding to what are presently the antifields of the antighosts, was 
discovered by Zinn-Justin in a discussion on the renormalisation of non-abelian gauge theories. 
It is therefore called the Zinn-Justin equation and can be considered as a rewriting of the 
Veltman-Ward identity for the BRST symmetry. 

4 In the sequel, in the quantum context, we will use simply 5* instead of S q to denote the full quantum action. 
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A detailed discussion of the BV formalism and renormalisation using the effective action has 
been given in [25|. Further advocacy for the use of BV in the formulation of renormalisability 
can be found in p6| . 



3 BPHZ 

To put the quantum theory of the BV method on a more solid footing, one wants a proper 
quantum version of the basic equations. In particular, because of its central importance, we 
are looking for a regularised version of the quantum master equation eq.(|2.6[) as an operator 
equation which is valid to all orders in perturbation theory. In this section we first discuss 
some alternatives, and then give our reasons for using the setup of Bogoliubov, Parasiuk, Hepp 
and Zimmerman. We continue with a summary of the basics of that method, and recall the 
treatment of the chiral anomaly paradigm. 



3.1 Alternatives 

As was already pointed out, the quantum modifications in eq.( |2.6D as compared to eq.(^ 
and in particular a non-zero AS", signal the possible appearance of anomalies. In fact, AS" is 
closely related to the Jacobian of the path integral measure under BRST transformations, with 



a potentially diverging determinant. Let us therefore consider the method of Fujikawa [27]. 
The BRST Jacobian determinant is a product of eigenvalues of some operator in an infinite 
dimensional space, and this is regularised by introducing a Gaussian damping, removing the 
contribution of the 'high-energy' modes. However, the method does not provide an unambiguous 
rule for choosing the regulator (the Gaussian damping), although the ad hoc conditions that 
are sometimes imposed to obtain a consistent anomaly can be replaced by a more systematic 
approach using Pauli-Villars methods j28|. Also, the extension of the method to higher loops is 
not straightforward. 

In contrast with the philosophy of the Fujikawa method, which considers the relevant Ja- 
cobian determinant as the primary object to regularise, a more comprehensive point of view is 
that a regularisation scheme should take care of the divergences of the full theory (up to a cer- 
tain order in perturbation theory). Most of the universally popular regularisation schemes are 
formulated as prescriptions to associate first a regularised, and in a second step a renormalised, 
finite expression to diverging loop diagrams. Let us then consider 'dimensional regularisation', 
which gives a method to replace the Feynman rules by an 'n-dimensionak generalisation, with 
n complex. Let us note that in this scheme one does not have at one's disposal a Lagrangian 
corresponding to the regularised step. It is well known that, when symmetry properties of the 
Lagrangian one considers change with dimension, dimensional regularisation requires a delicate 
treatment to provide the same answers as other regularisation methods, notably when comput- 
ing anomalies. In fact, in the initial BV-treatment dimensional regularisation was invoked to 
drop the anomalous terms. The lack of a Lagrangian seems to prohibit the formulation of the 
basic local BV equations in this scheme, although if one accepts the customary interplay of 
integer and complex dimensions one should refer to the treatment^] in M. 

We now turn to schemes that do have a straightforward Lagrangian interpretation. Originally 
due to Pauli and Villars, the simple device of subtracting from loops a compensating loop in 
which a particle with a very large mass circulates is the oldest among these. One writes down a 
Lagrangian for the original theory plus an additional piece where all terms involve one or more 
new massive (Pauli-Villars) fields. This is done in such a way that at this stage the theory is 

Especially charming through its use of d — 4 as an extra (global) variable of the theory, with a corresponding 
BV-antivariable that plays a role in the treatment of the anomaly too! 
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(one-loop) finite. Then one computes the effective action, and does the renormalisation. For a 
renormalisable theory this step absorbs the terms that would potentially diverge in the next step, 
where one lets all masses of the PV-fields tend to infinity. The contents of the original theory is 
defined to be what one obtains in this limit[| In this treatment, possible complications arise if 
the mass terms do not respect all symmetries of the original theory. Then the classical master 
equation is broken and one has to investigate carefully whether the breaking terms survive after 
taking the infinite mass limit. If they do, the theory may have an anomaly. This method goes 
together well with the BV formalism, and this way of computing the one loop anomaly is by 
now well-documented in the literature. It fits well because the method is based on a Lagrangian, 
just like the BV setupQ. 

The Pauli-Villars scheme only works well at one loop. One could imagine circumventing 
this restriction by using higher (covariant) derivatives in combination with the one loop PV 



scheme [29 ]. However, constructing such covariant derivatives requires an understanding of the 
geometry behind the gauge symmetry. For some models, for instance the W3 gravity model 
which we use as an example in this paper, such knowledge may be lacking. 

Another resolution to the divergence problem, operationally very close to Schwinger's proper 
time method, consists in modifying the propagators to cut off the loop momentum integrals. The 
Lagrangian gauge invariant implementation of this idea is necessarily nonlocal. This method 
has been advocated as a universal solution to the regularisation problems, and was applied 
to the present problem in ||. The most striking virtue of this method, namely regularisation 
to all loop levels, brings also the computation of higher loop corrections within reach JH], and 
makes it into the primary alternative to the method applied in this paper. We refer to (5]^[6| for 
a detailed exposition. 

3.2 Why BPHZ? 

In view of the coherence of the BV formalism, which rest very heavily on the use of a 
classical action functional, a properly defined quantum action functional is highly desirable, and 
preferably well defined functional integral expressions as well. When regularisation is formulated 
directly in terms of diagrams however, the path integral is just a shorthand notation for a 
collection of diagrams. This notation naturally suggests relations between correlation functions 
(=sets of diagrams), such as those obtained by functionally integrating by parts, the most 
familiar of these being the Schwinger-Dyson equations and Ward identities. Relations of this 
type are used heavily to set up the formal BV scheme, so we require that regularisation keeps 
them in manageable form. 

Of all renormalisation methods, the BPHZ method is probably the most solidly founded 
one. In fact it is used as a touchstone (see for example |32]) for other methods to establish 
their respectability. It adresses head-on the key issue in the regularisation of quantum field 
theory, viz. the definition of products of local operators at the same space-time point. In BPHZ 
this is done via an explicit prescription of all their matrix elements. These matrix elements are 
defined in terms of Feynman diagrams. The diagrams are the same as when working formally, by 
which we mean disregarding the divergence problem. When applying the Feynman rules as they 
follow from the action, one obtains meaningless expressions. The BPHZ method is a definite 
prescription, order by order in perturbation theory, to supplement the Feynman rules so that one 



6 The 5(0) problem of eq.(2.4) does not arise, because the offending right hand side can be avoided completely 
Jll| ] due to a cancellation between the original fields and the PV fields. 

7 We do not wish to distinguish here between Lagrangian and Hamiltonian methods, but rather set apart these 
methods from a purely diagrammatic approach. We will continue to use the term 'Lagrangian' to indicate both 
variants of action-based treatments. 

8 The prescriptions given in EjJ are not quite correct, but have been amended in the recent literature, see f3(i||. 
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obtains meaningful expressions. Whereas the prescription is inspired by formal manipulations, 
'subtracting infinite counterterms', in fact it makes no reference to such a dubious step. Let us 
remind the reader that the best-founded proofs of renormalizability of various theories rest on 
the application of this method. Let us also remind the reader that its prescription is universal, 
in that it applies irrespective of the order in perturbation theory. 

The big asset of the BPHZ scheme is that renormalized correlation functions with insertions 
of arbitrary composite operators can be defined in the same way and with the same ease as 
diagrams related to Green's functions of elementary fields. The technical tool that does the job 
is a so-called 'normal product' operator, that we discuss in the next subsection. The normal 
products uniquely determine the finite expression associated with every diagram. The relations 
between correlation functions containing composite operators that one formally derives from 
path integrals, all have counterparts as relations between normal ordered products, and there- 
fore as relations between renormalized correlation functions. Thus while the BPHZ method is 
not exclusively based on an action functional, it does provide well-defined quantum operators 
that correspond to classical functionals through the normal product concept, and it does pro- 
vide counterparts to the classical equations obeyed by these functionals. The BPHZ scheme is 
therefore well-suited to formulate the basic local equations of the BV formalism. 



3.3 Basics of BPHZ renormalization 

Elaborating on the work of Bogoliubov, Parasiuk and Hepp [33|, Zimmerman 34] set up a 
renormalization scheme in which renormalized Green functions with composite operator inser- 
tions are obtained in essentially the same way as renormalized Green functions of elementary 
fields. The method^, which bypasses the regularization step, determines the renormalized coun- 
terpart Ir{p) of a divergent loop integral I{p) associated with a specific one particle irreducible 
Green function by subtracting from its integrand the first few terms in the Taylor series around 
vanishing external momenta p, i.e. 

T(p) = Jdk l(p, k) — ► I R (p) = J dk (l - t 5 p ) T(p, k) . (3.1) 

The minimal value of 5 is the superficial degree of divergence of the diagram (if the diagram is 



not finite by power counting). One- loop diagrams can unambiguosly be treated with eq.(|3.1|), 
but multiloop diagrams with overlapping divergences require the use of the forest formula p3] 
to determine the Taylor series to be subtracted (see section ^ for an explicit example). 

In fact, the BPHZ renormalization scheme allows to define for any local composite operator 
O(x) a sequence of normal ordered composite operators 

N a [0{X)\ , 

with a an integer greater than or equal to the canonical dimension of the operator 0{x) (counting 
1 for every derivative, etc. ). A normal product is called minimally subtracted (resp. oversub- 
tracted) when a is equal to (resp. greater than) its canonical dimension. For a generic connected 
Green function of free field normal products N Qi [Oj(xj)] 

^Y[^a l [O l (x l )]^ , 

the order 5 of the Taylor series in Zimmerman's prescription ( |3.1| ) is given in terms of the degrees 
a, as 

S = n + ^( Gi - n ) - , (3.2) 



3 Useful general references, containing a more detailed exposition are [^, |ll| . 



8 



where n is the space-time dimension and the last term subtracts all the canonical dimensions 
dk of external lines. 

The use of oversubtracted operators, while seemingly superfluous, is in fact one of the virtues 
of the BPHZ method and a cornerstone in its use for anomaly computations. Moreover, it 
is useful to introduce not one but a variety of oversubtracted normal products p^] . Indeed, 
the degree of oversubtraction of an operator can be propagated to subdiagrams^j in different 
ways, the only constraint being that the degree of oversubtraction of a subdiagram should be 
at least as large as the degree of oversubtraction of its subdiagrams. When subdiagrams are 
always maximally oversubtracted one speaks of isotropic normal products (all lines are treated 
in the same way), otherwise of anisotropic normal products. We do not need the most general 
anisotropic ones, but only the Gomes-Lowenstein pjl (GL) version, 



N dl+d2+Q [{V X }V 2 ] , (3.3) 
for which the rule that determines the subtraction degrees 5 7 of subdiagrams 7 is modified tcj^| 

( di + d 2 - n, if i G 7 and L(V 2 ) C E(j) 
5 7 = n -^2d k + ^2(ai - n) + < d\ + d 2 + a-n, if i £ 7 and L(V 2 ) ^(7) > (3-4) 
k i^io I 0, if io ^7 

where io is the normal product vertex under consideration; di the canonical dimension of the 
operator Vi\ L(V 2 ), the set of lines associated with V 2 and £"(7) the set of external lines of 
the subdiagram 7. Note that this rule implies that if V\ is linear in the fields then there is no 
distinction between the isotropic and the anisotropic products. 

Let us conclude this summary by indicating two fundamental and useful properties of the 
normal ordered composite operators. The first of these, the BPHZ renormalized version of the 
Schwinger-Dyson equations is0 

m / N a p{x)\ • • • ® Am i 



5<S> B (x) 

N. 




' a+n- 



5<S> B {x) 1 K n 5<Z> B {x) 



^( Xl )...^(x m ) , (3.5) 



where is the engineering dimension of the field $ B (x); So, the massless free part of the action 
and Si lumps together eventual mass terms and interactions, i.e. Si = Si + S m . Notice the 
anisotropic normal ordering in the second term on the rhs. 

The second property, the so-called Zimmerman identity, establishes linear relations among 
different normal products of a given monomial. Different formulations may be given. A first 
version relates normal products with different subtraction degrees a\ > a 2 to all operators Oi{x) 
having the same quantum numbers and with canonical dimensions < ct\ 

N ai [0(x)\ - N Q2 [0{x)\ =J2 r * N «i Pi(x)] , (3.6) 



I We always refer to proper subdiagrams, i.e. subdiagrams that are themselves 1PI. 

II The Taylor series operator acts on all momentum factors produced by the Oi, but not on the propagators of 
the lines that are external to the subdiagram 7. 

12 Here and in the following, to avoid unnecessarily heavy notations, we act as if there are no dimensionful 
parameters in the model, and all operators in the action are minimally subtracted N n [ ■ ]-products. It is rather 
trivial to overcome this restriction. 
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where the coefficients rj vanish at lowest order. Our approach in fact uses a consequence of 
( |3.6[) , relating oversubtracted operators to their minimally subtracted counterpart 

N a+Q [0(x)\ - N a [0{x)\ = N ai [Oi{x)\ , (3.7) 

i 

where now all the operators on the rhs are minimally subtracted. A second, extended version 
of the Zimmerman identity further relates anisotropic and isotropic normal products with the 
same substraction degree 

N dl+d2+a [{Vi}V 2 \ - N dl+d2+a [Vi Vi\ = r ^+d2 [Oi{x)V 2 ] , (3.8) 

i 

where the same restrictions for and apply also in this case. In the next section, the two 
relations ( |3.5|) and ([T?]) will be used extensively for the implementation of BPHZ ideas into the 
BV framework. 



3.4 Anomalies and the massless case 



The exposure of the chiral anomaly in the present framework is an application of the BPHZ 
method that dates back some 25 years. We shall not repeat this calculation here, since a short 
account of it can be found in ||, but just sketch the main ideas and stress that such recipe 
works in general for determining the anomaly in a global, continuous symmetry. 

The central equation, which computes the divergence of (matrix elements of) the axial vector 
current in Q.E.D. (@|, eqs. (1.7) and (1.8)) is 



a„ ( n 3 



(75^(a0) 



5^(x) 



ijyf^ij} (x)X 
5 



5ip(x) 



X ) + 2im ( N 4 V>7 ^ (x)X 



(3.9) 



which is an application of eq.( |3.5| ). Then the oversubtracted operator N4 [i/>7 5 'i/>] is rewritten in 
terms of minimally subtracted ones by using Zimmerman's identity ( |3.7| ) as 



Ni 



N, 



+ sN 4 



(3.10) 



An important feature of the treatment is that, for massless fermions, a mass has to be 
introduced to perform the calculation. The reason is that the BPHZ method introduces an 
apparent infrared divergence in the subtraction terms, in spite of the fact that no physical 
infrared problems are associated with the vanishing of the electron mass. It comes up when one 
tries to calculate the matrix elements of the oversubtracted operator N4 [t/>7 5, 0] : they diverge 
in the massless limit. Ultimately, this is why one is forced to reduce the oversubtracted to an 
ordinary normal product: the latter has matrix elements that stay finite when the mass goes to 
zero (if there is no physical infrared problem). The upshot in the present computation is that, 
in eq.( |3.9| ), when m — > 0, one picks up the terms from eq.( |3.10| ) that go as m~ 2 . This is the 
origin of the axial anomaly. For the actual computation we refer to the quoted literature for the 
axial anomaly, and for other examples to sections 5.3 and and appendix of this paper. 



4 Anomalies in the BV— BPHZ renormalized framework 

The properties of the normal ordered products (^), ( |3lj| ) have been known to provide 
a natural framework for the investigation of anomalies in global symmetries for a long time. 
Furthermore, they were extensively used in the original presentation of the BRS quantization 
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scheme [16]. We will carry these ideas one step further by implementing BPHZ ideas in the 
BV framework, in order to formulate its basic local equations -essentially the quantum master 
equation- on the level of renormalised perturbation theory and to derive a proper method to 
compute higher loop anomalies. For treatments on the level of the effective action, which have 
been limited to one-loop anomalies, we refer to []ll| and references therein. 

Roughly speaking, this program is realized by substituting formal products of fields by 
suitable Zimmerman normal products. In the assignment of dimensions d A , d* A to the fields 
<& A and antifields & A , there is a certain degree of arbitrariness, although they are largely fixed 
by the classical free Lagrangian. This arbitrariness is fixed in a way that is convenient for the 
BV framework by demanding that the BRST operator preserves canonical dimensions^. The 
following relations then hold 



d 



SS 



5<5> A 



n 



d: 



d 



5S 



df 



n 



d 



A • 



(4.1) 



such that the dimensions of a field and of its corresponding antifield add up to the space-time 
dimension n. As a consequence, normal products involving BRST sources obey, for example 



N a 



A 1N a— n+d.A 



R 



(4.2) 



An essential ingredient to deal with massless theories within our program consists in the 
introduction of a small mass m for the massless fields. As was pointed out in sec.^J, this is 
necessary to prevent the appearance of spurious IR divergences when performing the prescribed 
subtractions at vanishing external moment aP|. In these cases, one therefore replaces the original 
massless action S by a massive action S 



5($,$*) 



m 



S($,$*) + ^S m ($) 



(4.3) 



supplemented with the rule that the limit m — > is taken at the very end of the renormalized 
computations. Also, d A + d B + d(TAB) = n — 2. 

We will now formulate the basic equations governing the BV method at the renormalised 
perturbative level. As already pointed out in the introduction, in the BV scheme the quantum 
BRST structure and its possible violations are naturally described by means of the BRST- 
Veltman-Ward identity (|2.8| ) for the effective action T. To write down the corresponding local 
equation we start from the generating functional of all connected Green functions W(J), of 
which r is the Legendre transform, r(<£, $*) = W(J, <&*) — J<3?, where it is understood that to 
obtain the left hand side one should solve the equation $ = ^j- for J as a function of $ and 
<&*. By inverting the Legendre transform, one then also has that = ~~ J- The left hand side 
of equation fl2.8|) is then equal to 



d n x J A {x) 



5T 



(4.4) 



13 If one follows the traditional assignment, of |9[ for instance, the usual Yang-Mills ghosts and antighosts would 
both have (scaling) dimension 1, whereas we assign dimension to c and 2 to 6, as in ref. Also, antifields 
would have dimension (being external fields), whereas our antifields do not. The BRST operator increases that 
(traditional) scaling dimension by 1 unit. We will adopt conventions matching those of conformal field theory 
(see table [l] on page [H] and table ^ on page ^J). The difference is only one of bookkeeping. 

14 Subtraction at non-zero momentum is feasable, but unnecessarily complicated regarding counterterms and 
Lorentz invariance. 
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A general property of Legendre transformations is that the derivatives of T and W with respect 
to any parameter, $* for example, are equal: 



<5r($,$*) _ 5W(J,$* 
5$* A (x) ~ 8<Z>* A (x) 



(4.5) 



where on the l.h.s is kept constant in taking the derivative, and on the r.h.s. J is kept constant, 
and it is understood that after taking the derivatives the functional relation between J and $ is 
imposed. Since W is nothing but the generating functional of all the connected Green functions, 
a (formal) power series in the source J, its derivative with respect to the parameters 3>* is the 
generating functional of all connected Green functions with one extra insertion of the operator 

~5s 



. Suppressing the $* dependence in the notation, we will denote the resulting generating 



functional with the insertion of X by (X) cJ , so that, for example, W(J, <£*) = (l) c j. This 
local functional of the fields (and possibly of antifields as well) that one has to insert is not 
just the classical functional (which would be ill-defined), but one of the normal products that 
correspond to this classical functional, which have a definite meaning in the BPHZ scheme. The 
particular normal product one takes is an integral part of the definition of the quantum action 
functional one starts from. We take the interaction terms in S, in particular the antifield terms, 
to be minimally subtracted. Other choices differ in terms of order Ti or higher, see (|3.7|). One 
may wish to include in the action more general terms of higher order in fi. This would be the 
case when loop anomalies are to be canceled by finite non-invariant counterterms (usually called 
Mi,M2, ... in quantum BV formalism). We will continue to denote the full quantum action, 
including these terms if they are introduced, by the symbol S. (If some of these terms contain 
antifields, they are to be viewed as quantum modifications of the BRST transformation laws.) 
The normal ordering to be used for the lhs is then, according to ( |4.1| ), N d A [■]. The equation 
(|2.8|) is then completely equivalent with 



d n x J A (x) N 



d A 



5S 



ih / d n x (N„ [A(x)]) 



C.J 



(4.6) 



C.J 



It should be emphasized that this equation is an exact one, involving on the l.h.s. only quantities 
which are well-defined in renormalised perturbation theory, although rather empty until one 
specifies the operator N n [«4(x)] on the r.h.s.: in fact one may take it as the definition of A. 
If one treats the equation formally, ignoring the normal product rules, for example using the 
heuristic path integral formulas, then, as already mentioned in section |2,3| , it is rather easy to 
derive the general formal expression A = AS + -^(S, S), for which there is of course no place 
in the BPHZ scheme. Reverting to a more proper treatment, there is still a general statement 
that can be made: a celebrated result of general BPHZ renormalisation theory, namely Lam's 
theorem ]22|, 23[ (also called the quantum action principle) guarantees that the anomaly insertion 
A(x) in ( |2.8 ) is the integral of a local composite operator, whose canonical dimension is n 
in our conventions 15 , without specifying its explicit form. We will not need to rely on this 
general theorem here, since we will now deduce an explicit formula for A that shows this locality 
explicitly. 

Expanding (4.6) in the sources J a one obtains 



d n x 



5S 



6<f> A (a 



-(-1) €A <S> Al ( Xl )...<S> y 



The well-known 75 anomaly in 4 dimensions, including the factor c, would have canonical dimension 5 if we 
had adopted the counting with the usual scaling dimensions. 
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Jd n x (N n [^(x)]d> Al (^i)---^ Am (^ m ) 



(4.7) 



These equations are nothing but the typical anomalous Ward identities for the usual BV-BRST 
transformation 5 since, for antifield independent functionals, one has classically that 



d n x 



5S 



6$* A (x) 5^ A (x 



i-iy A \f[$] = (f,s). 



(4.8) 



To obtain from eq.( }4.7D an equation in terms of local quantum operators (instead of corre- 
lation functions), we apply the BPHZ expression of the Schwinger-Dyson (SD) equation, in the 



5s 



This 



form given by Gomez and Lowenstein, ( |3.5| ), to its lhs, for the operators O(x) 
leads to the following operator identity in terms of (possibly oversubtracted) anisotropic normal 
products 



-ih J d n x N n [A(x)\ 



d n x Nr. 



5 S 5S 
5<S> A (x) 5<f>* A (x) 



+ N n 



<5(S m + S 7 ) SS 



5$ A (x) I 5<f>* A (x) 



Assuming the validity of the classical master equation (2.3) for S = So + Si, the functionals 
inside the normal ordering can be simplified, resulting in 



-ih J d n x N n [A(x)\ 



d n x N n 



5S 



6_Si_ 
5<f> A (x) | 5$* A {x) 

5 S m 



5S! 5S 
5§ A (x) 5<$>* A {x) 



+ d n x N n 



5S 



5$ A {x) I 5<5>* A {x) 



(4.9) 



This is the promised explicit local expression for the BRST anomaly in the BPHZ frame- 



work. We now comment on the different contributions in (4.9). The first contribution, poten- 
tially present in any interacting theory, comes from the difference between the anisotropic and 
the isotropic normal products of the BRST variation of the interaction term Si, as a straight- 
forward application of the generalized Zimmerman identity ( |3.8D shows. These contributions 
vanish however for minimally subtracted^ BRST variations of the interaction, since then the 
GL anisotropic normal products (|3.3| , |3.4|) and the isotropic ones exactly coincide. The second 
contribution to the anomaly originates in the explicit breaking of the BRST (or gauge) symmetry 



by the IR regulating mass term S m in (4.3) 



lim m 2 I d n x N r , 

m->0 



® a Tab 



5S 



(x) 



(4.10) 



The massless limit may result in a non-zero contribution, as can be seen through a generalization 
of the argument presented in sec. 3.4. The composite operator in ( 4.1C| ) has dimension n — 2 



16 In [|35[ it is indicated that in nonlinear cr-models one may prefer non-minimal interaction terms. Therefore, 
these terms can not be dropped in general, although in the examples of this paper they could. See also footnote [12 
on page H 
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(or n — 1 for Fermi fields, for which the typical mass term in (^^) would be proportional to m 
instead of m 2 ). Due to the normal ordering degree n, its insertion in correlation functions leads 
to over subtracted integrals, so that taking the limit m — > naively by simply discarding such 
contributions is incorrect. To compute this limit correctly requires a conversion to minimally 
subtracted operators first. We do this conversion in two steps. First we use the Zimmerman 



identity (3JS) to relate the anisotropic, oversubtracted normal product to its isotropic, oversub- 
tracted counterpart, where the normal product corrections are expressed in terms of operators 
Aik(x) with canonical dimension dk < n, i.e. 



$ A Tab 



6S 



(x) 



B 



5$ 



B 



+ £^N dfc [M k (x)) . 
k 



In a second step, the isotropic, oversubtracted operator is expressed in terms of minimally 
subtracted normal products by means of the second Zimmerman identity Q3.7| ) . The net result 
of this process is then written as 



N, 



$ A T AB 



SS 
6$* B 



(x) 



+ $> fc N dfc [M k (x)\ 



(4.11) 



where the numerical coefficients are at least of first order both in coupling constants and in 
ti. After substitution of ( |4.11 ) in ( 4.10| ), the limit m —* can be taken. If the theory has no 
infrared singularities the minimally subtracted normal products have finite matrix elements and 



consequently only those terms will contribute to eq.( 4.10 ) that have coefficients b^ of order 1/ 



Ira 



The corresponding operators determine the second contribution to the anomaly in eq. fl4.9j ). 
If there are infrared singularities, the effective action T strictly speaking does not exist. A 
treatment of a Wilson-type effective action can certainly be developed, and a corresponding 
local operator treatment, and it seems reasonable to anticipate that, as long as the infrared 
divergences are merely logarithmic the treatment given above remains valid 

To conclude our theoretical approach, we remark that the contribution to the anomaly 
generated by the IR regulating mass term is not uniquely fixed. Indeed, there may be more 
than one way to construct a mass term in ( |4.3| ), Different mass terms may behave differently 
under some symmetries, causing the anomaly to show up in different symmetries. A similar 
ambiguity is present in the context of the one- loop Pauli-Villars regularisation (see [28, fl, and 
[j36f for an example how this can be exploited), but we have not investigated whether this 
correspondence is complete in all details. 

In summary, the simple rule of substituting formal products of fields by suitable Zimmerman 
normal products has led us to an explicit expression, eq. ([4.9D , for the anomaly in the BV-BPHZ 
renormalized framework, which is theoretically quite interesting but not necessarily practical. 
In fact, it turns out that for the actual computation of the BRST anomaly insertion A in (|2.8j ), 
eq.([Oj) is more convenient once an enumeration of possible monomials in A has been given. 
The next section will illustrate this by computing the conformal anomaly in chiral W% gravity 
using both alternatives. This will of course also serve as an illustration for the more abstract 
development of this section. In the section after that, we will go on and apply the method to 
the one and two-loop anomaly of W% gravity, mainly using the anomalous Ward identities 



but the two-loop anomaly will also be derived from eq.(4.9) 



5 The anomaly of chiral W2 gravity 
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5.1 The model 



Chiral W-Vgravity is a 2-dimensional model (n = 2) describing D matter fields cp 1 , i = 
1,...,D, coupled to a gravitational gauge field through their energy- momentum tensor T = 
l/2(dcf) t )(d(f) t ). A convenient, gauge- fixed extended actiorf 7 ! for this model can be taken from 

m : 



S 



I 



d 2 x 



1 



(d0 l )(cV) + &(dc) 



+4>* [c(50*)] + b* [-T + 2b(dc) + (db)c] + c* [(dc)c}} 
S + J d 2 x <5>* A (x)R A (x) . 



(5.1) 



In writing (5.1) the following conventions have been used: 



d = d 4 



B = d. 



V2 



(x 1 ± x c 



The fields c and b in ( |5.1| ) are respectively the ghost and antighost of the spin j = 2 gauge 
symmetry, and {(f)*, c*, b*} are the antifields or sources for the BRST transformations for all fields. 
Remarkably, in this formulation interactions are completely contained in the antifield dependent 
terms, so that BRST sources can be regarded as coupling constants, in which perturbative 
expansions can be performed. In fact, the antifield b* appears in (|5.l| ) as the source for (minus) 
the total energy-momentum tensor of matter and ghosts, viz. 

— (dPXdfi) + 2bdc + (db)c = —T — r gh = -T t . (5.2) 



R b 



Thus b* takes over the role of (minus) the gravitational gauge field. 

We will work with the more general expression for the ghost energy momentum tensor 

Tj = -jbdc+(l-j)(db)c, (5.3) 

for a generic spin (j, 1 — j) 6c-ghost system. At the moment we need j = 2; in the next section, 
for chiral W3, we will have both spin j = 2 and spin j = 3 ghost sectors. 

For the assignment of dimensions, we follow common practice in conformal field theory. There 
are two conformal dimensions for each field, corresponding to the left and right Virasoro algebras. 
The total dimension d is the sum of these, the spin j is the difference. The relevant assignments 
are collected in table |l[ For the antifields, we extended the assignments in accordance with 
eq.([4.1|), and spin conservation. Note that the fields and the d derivative have vanishing d — j, 
whereas antifields and the 3 derivative all have d — j = 2. This will be useful in our analysis. 

These dimensions do not correspond to the assignments used in section [D| a canonical 
assignment would attribute dimension 1/2 to all ghost fields. The subtraction rules in eq.( |3.2| ) 
and eq.(|3,4|) require modifications that reflect this. We quote the resulting modification only 
for the case of operator products that do not contain antifields (which should be extracted first 
using eq.(O)), both for W2 and W3: 



5 = n + ^2(ai - n) + ^ (4 - 1) , 



(5.4) 



where the primed sum runs only over ghost external lines. This new rule extends to the 
anisotropic subtraction degree (3^), in which J2k^k should be replaced by J2'k(~ dk + 1) as 
well. 



17 We take this at the same time to be the quantum action, and as in the general treatment of section ^ the 
minimally subtracted normal products are understood. With this understanding, and including no modifications 
Mi of higher order in fi, it will turn out that the anomaly comes out in its conventional form. 
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dim d 


spin j 


gh-tt 


d-j 
















c 


-1 


-1 


1 





b 


2 


2 


-1 







2 





-1 


2 


c* 


3 


1 


-2 


2 


b* 





-2 





2 


d 


1 


1 








8 


1 


-1 





2 



Table 1: Additive charges of fields and antifields in the W2 model. 



Finally, for the reason mentioned in section 3.4, we introduce a mass term for all propagating 



fields. For the matter fields, this is conveniently done by including in S m (|4.3[) the term 



2 

771 



- J d 2 x (f) i {x)4) i {x). (5.5) 
The free propagator for the massive matter fields then becomes 

( ) = A«(x -y) = & * h 2 8(x - y), (5.6) 

\ / 00 — m z 

with the subscript indicating from now on free correlation functions, that is, the ones evaluated 
in the zero-th order in the perturbative expansion in antifields. Introducing a mass term for 
the ghosts is slightly more subtle. The obvious choice, — m J d 2 x be, has the inconvenience of 
breaking spin (j) invariance. This can be amended by using the non-local mass term 

— m 2 [ d 2 x b(x) —c(x) . (5-7) 

J a 

The introduction of this nonlocality requires some comments. The extra mass terms recall a 
similar feature in Pauli-Villars regularisation. There, extra fields are introduced to regularise 
loop integrals, and these extra fields are given very large masses at the end of the calculation. 
Anomalies are understood as remnants of the fact that these mass terms do not respect all 



symmetries that were present classically. For example, in [28] the ghost number anomaly was 
traced to a mass term with ghost number different from 0. There are some important differences 
with the BPHZ method followed in this paper. The most important one is that now the masses 
are given to the original fields (there are no extra fields), and they tend to zero at the end, 
not infinity. Thus the nonlocality of the mass term does not pose any problem in the present 
framework, whereas in the PV scheme it would. The locality is preserved here by the subtraction 
method (a Taylor series in momenta), and the introduction of masses is a device needed when 
oversubtracted operators appear. The consistency of the setup can be verified by checking that, 
despite having ghost number zero, this mass term reproduces correctly the well-known ghost 
number anomaly. This is demonstrated in appendix [C]. In fact the peculiar non-locality is at the 
origin of the ghost number anomaly. From a perturbative point of view, the nonlocal character 
of (|5.7| ) poses no problem either. With this choice, the free propagator for the massive ghost 
fields 



(c(x)b(y)) = (b(x)c(y)) = - (b(y)c(x)) 



G(x-y) = -jB^5(x-y). (5.8) 
00 — m z 
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takes a form that is common for Fermi fields. We conclude that the ghost mass term (|5.7|) , 
though non-local, is indeed appropriate for our purposes. 



5.2 The W 2 anomaly from the anomalous Ward identities 



Now everything is ready to start the computation of the anomaly by using the anomalous 
Ward identities (4/7). Its evaluation can be divided in two steps. First, one uses all possi- 



ble information concerning symmetries and quantum numbers to determine the general form 
of the monomials that can build up the anomaly. The connected character of the involved 
Green functions together with the presence of the normal products severely restricts the form 
of the anomaly candidates, so that their field content is completely determined. A standard 
perturbative computation then determines the coefficients of such anomaly candidates. 

Let us now implement this program in the present example. Since the BRST transformation 
5, eq,(fO|), preserves spin (j) the anomaly is a spin zero object, j[^4(x)] = 0. Combining 
this information with the canonical dimension, d[.A(x)] = 2 in our conventions, one finds (d — 
j)[A(x)] = 2. From Table |l| it follows that all possible terms in the anomaly contain either 
precisely one antifield or precisely one 9-derivative. Lam's theorem further ensures that the 
numerical coefficients of the candidate anomalies are at least of first order in both U and the 
coupling constants, the antifields <£*. Thus the anomaly has no d dependent terms, and is of 
the form A{x) = &\(x)F A (d, x). This determines the antifield dependence of the anomaly 
completely. 

Consider now the anomalous Ward identities ( ^4.7[ ). From the the above analysis, it is clear 
that all the relevant information about the functions F A (x) -up to now, still arbitrary functions 
of d and all fields- is contained in the loop contributions to the linear term in the antifield 
expansion of the lhs of (|4.?1) , namely 



d 2 y( N dA R A (x) \N dB R B (y) 



ih( N 



d A 



5$ B (y) 
F A (x)] <& Al {xi)...<S> Am {x 



iy B <$> A i( Xl )...<$> Am (x r , 



c.O 



c.O 



where relation 



(5.9) 

can now be completely deter- 



was used. The form of the functions F A (x 
mined as follows. On the one hand, the normal ordered character of the insertions in (5.9) and 
the connectedness of the Green functions indicate that "proper" loops (i.e. no tadpoles, which 
are zero due to the subtractions) can only be generated by making at least two contractions 
between the two BRST transformations R A and R B on the lhs of (^9|)0. Their quadratic field 
dependence restricts moreover the maximum number of contractions to be precisely two and 
leads thus to the conclusion that only one-loop anomalies can appear. On the other hand, from 
inspection of the free propagators (|5.6|) , (|5.8|) , it follows that these double contractions are only 
possible when both A and B refer to the antighost b. These contractions made, no fields are 
left to contract with, such that the "test product" of fields & Al . . . <& Am must be equal to b(z) 
in order to get a non-vanishing lhs. The lhs of fl5.9j) becomes 



f d 2 z ^ N 2 


R\x) 


N 2 


R b (z) 











since R b 



-T - U 



-T t , see eq.(|5. 



5b(z] 



(N 2 [T t (x)]N 2 [T t {y)\ 



c,0 ' 



c.O 



The correlation function on the rhs of (|5.9| ), of the form ( N 2 [-F(x)] b{z) ) c with the above test 
product, is non-vanishing if and only if F{x) contains precisely one c-ghost, since the presence of 

18 As a consequence linear transformations do not contribute. 
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Figure 1: The one- loop anomaly diagram for W 2 . 



other fields in F(x) would lead to extra tadpoles. Spin conservation implies A(x) 
ab*d 3 c, so that (5.9) becomes 



ab*F{x) 



<N 2 [T t (x)]N 2 [T t {y)\ 



c,0 



iha{N 2 d 3 c(x) %)}c,0- 



(5.10) 



This is an explicit equation for the anomaly coefficient a, which corresponds diagramatically 
to Figure [I]. The coefficient a is thus determined by double contractions of the total energy- 
momentum tensor T t with itself. Let us stress that the lhs of ( 5.10 ) is completely determined in 
renormalized perturbation theory, owing to the normal orderings. The equation expresses in the 
BPHZ language the well-known fact that double contractions in the OPE of Tt(x)Tt(y) generate 
the conformal anomaly. 

The anomaly computation is rounded off with the determination of the actual value of the 
coefficient a from eq. ( |5.10| ). On the lhs, the matter and ghost field contributions can be analyzed 
separately, since no contractions between T and T g h are possible. The formal expression for the 
contribution of the matter sector is 



(T(x)T(y)) 



,D 



2 



dd — m 



■.6(x - y) 



-Ti 



2 D f d 2 p 



where I 2 2(p, m ) is a quadratically divergent integral (the explicit expression is in eq. ( gig ) with 
a = b = 2). The correct expression is similar, and arises by replacing the composite operators T 



by their (minimally) normal ordered counterparts, N 2 [T]. The subtraction rule (5^) states that 
the integrand 222(1, P, m ) i n t ne quadratically divergent one-loop integral must be replaced by 
(1 — t 2 ) 122(1, P, m ), confirming that minimally subtracted operators yield subtraction degrees 
precisely equal to the superficial degree of divergence. The correlation function of the normal 
ordered energy-momentum tensor for the matter sector is 



<N 2 [T(x)]N 2 [T(y)]) cfi 



-Ti 2 



D f d 2 p 



(2ny 



Jp-(-y)^-t 2 p )-I 22 (p,m) 



where we introduced the short hand notation 
lg(p,m) = (l-t n v )-I ab (p,m) 



J^k^-^^P, 



(27T) 



and I a b(q,p,m) is the integrand of the formal integral 
d 2 q <f(p — Q) h 



Iab(p, m) 



(2ir) 2 (qq + m 2 )[(p - q)(p - q) + m 2 ] 



d 2 q 

(2TT) 2 



l a b(q,P, m) 



(5.11) 



(5.12) 



(5.13) 
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with a, b > 1. For general spin j the formal expression 

(2j 2 - 2j + l)/ 22 (p, m) + (2j 2 - 2j)I 31 (p, m) 



^ P 4p-(x-y) 



(T 3 (x)T J (y)) cfi = h 2 j j±* 
is likewise converted to the well-defined 



d 2 p 
(2tt) 2 



Jp-(x-y) 



<N 2 [T^^Na [T j (y)]) c>0 = 
(2j 2 - 2j + 1) (1 - t 2 ) • / 22 (p, m) + (2j 2 - 2j) (1 - t 2 p ) • / 31 (p, m)] (5.14) 



The explicit computation of the BPHZ renormalized expression I$,(p) is straightforward, 
and details are relegated to appendix^. We simply quote the final result in the limit m —> 0, 
namely ( |A.9| ) 



,o+6-l 



i2 = ~ ^-B(a,b) 
Ztt p 



(5.15) 

with B(a,b) the Euler Beta function. The renormalized correlation functions corresponding to 
(ETllD , ( gig ) are found to be 



<N 2 [T(x)]N 2 [T(y)]> c , = 



2 ,. j2 

24vr 7 (2tt) 2 
i?i 2 r 



^ 2 -P f d 2 ^ c ip.{ X -y)lP 
P 



2(6j 2 - 6j + 1) 



d2 P i P -(x-v)P 3 
(2tt) 2 p 



(5.16) 



(2tt) 2 p 



and the renormalized expression for the lhs of ( 5.10 ) is 

<N 2 [T,(x)] N 2 [T t (y)]) c>0 = ^-(D - 26) J 
After computing the rhs of ( |5,10| ) in the momentum representation in the limit m — > 

/d 2 p 
(27T) 

a comparison of ( 5.17] ) and ( |5.18| ) gives 



(5.17) 



d2 P c ip-(x-y)P 3 
p 



(5.18) 



a =2T^D-2,). 

Combining this with the previous results about the functional form of the anomaly, one finally 
obtains 



A 



i(D - 26) 
24^ 



d 2 x cd 3 b* , 



(5.19) 



a well-known result. 



5.3 The W2 anomaly from Zimmerman identities 

In the previous subsection we showed the relation of the anomaly with the Ward identity for 
correlation functions. In this subsection we show its relation with the local operator equation 
( [4 . 9| ) . It may arise from two different sources: normal ordering corrections originating in the 
anisotropy of the variation of the interactions, and the symmetry breaking produced by the IR 



19 



regulating mass term. If we take both the interaction terms and the transformations in ( |5.1[ ) to 
be minimally subtracted normal products, only the second source contributes: 



lim 



d 2 x 



No 




Ax)} 



N 



SSr, 



5S 



S<f> A (x) 5<f>* A (x) 



(5.20) 



We have inserted the term with the minimally subtracted operator, which is supposed to vanish 
in the limit m — > anyway, to bring out the fact that classically this quantity is zero. The 
variation of the mass terms (|5.5|) and (|5.7|) reads explicitly 



N 2 



5S 



5<S> A (x) I 5<f>* A (x 
+m 2 N 2 



-m 2 N 2 



^{c^} (*) + (-<:) { / }( ) ) 



1 



{{dc)c\ (x) + Qc ) {2bi)r + U)b)r\{:r) 



(5.21) 



A non-zero result for eq.( 5.20| ) arises from a difference in the subtractions. We first spell out 
the influence of the anisotropy on the various terms in eq.( 5.21 ), and later turn to the influence 
of the oversubtractions. For the terms involving matter fields, when inserted in loop correlation 
functions^, the absence of interaction terms containing the product bft forces the factor ft to 
be always connected to an internal line, whereas the factor (^cj , instead, is always connected to 
an external line. Both then produce matter loops only. As a consequence of the definition ( ]3.4| ). 
one sees that the first of these anisotropic normal products behaves effectively as an isotropic, 
oversubtracted one, whereas the second reduces in fact to its isotropic, minimally subtracted 
version (and thus produces no contribution to the anomaly (|5.20Q ). Similar conclusions can be 
drawn for the first of the pure ghost terms of eq.(5.21). Since the factor (hb) is necessarily 
internal in a loop, the anisotropic normal ordering effectively behaves once again as an isotropic 
oversubtracted product. For the second pure ghost term, the net effect of the anisotropy is to 
first eliminate ghost loop contributions formed by considering the factor (hcj as external, and 
to force the remaining ones to be isotropic, although oversubtracted. 

Having disentangled the consequences of the anisotropy, the oversubtraction can now be 



taken into account. This is done through the Zimmerman identity (3.7) 



N 2 [9(x)] - N [6(x)] 



Y,Pk^d k [M k (x)], 

k 



(5.22) 



with d[A^fc] = dk < 2, and with the composite operator 6(x) given by 



@(x) = e m (x) + G gh (3 



-ft\cdft{x)\ 



+ 



1 



g bj {{dc)c} (x) + (JcJ {2bdc+(db)c} (3 



(5.23) 



The coefficients are most easily determined by inserting both sides of ( |5.22 ) into correlation 
functions with arbitrary products of fields. The insertion 0(x) (|5.23| ) has spin and dimension 
zero, which restricts the normal product corrections on the rhs of ( |5.22 ) to be exactly of the form 



m 



with f>k mass and coupling constant independent numerical coefficients. 

and elementary fields is 



19 The only non-zero correlation functions that can be formed with ( |5.2l| 
(N [A] (x)b(y) ) Q . See the previous subsection. 
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The actual value of the coefficients pk and the field content of the objects F A is again contained 
in the loop contributions to the term linear in antifields, namely 

d 2 y / N 2 [<Z> A (x)R A (x)] {N 2 [@(y)] - N My)]} ]J ^(x { ) ) = 

5>* I d\ /n 2 h* A (x)F k A ( X )} n^(^) ) • 



n idx 



c,0 



(5.24) 



c.O 



Considering possible double contractions between Q(y) and R A (x) = R b (x), which always leave 
a factor c(y) free, the test product of fields is again b(z) and the operator on the rhs can only 
be <&* A (x)F A (x) = b*d 3 c(x). Differentiating w.r.t. the antifield b* , ( |5.24 ) reduces to 



d 2 y ( N 2 [T t (x)\ {N 2 [@(y)} - N [@(y)}} b(z) ) cfi 
p(N 2 [d 3 c(x)] b(z)) c0 , 



(5.25) 



which is again an equation for the numerical coefficient p of the potential anomaly A = b*d 3 c. 

Let us now evaluate the contribution p m of the matter sector -coming from double contrac- 
tions between Q m (y) ( |5.23| ) and T{x)~ to the coefficient p in eq. ( |5.25| ), The BPHZ renormalized 
correlation function is 



K(x -y) = 



d 2 z (N 2 [T(x)]N a ^(z)c(z)(^)(z) b(y) 



c,0 



h 2 D 



d2 P Jp-(x-y) 



(2nf 



pp + m- 



:(l-t;)-I n (p,m) 



The short hand notation ( 5.12 ) was used again. The property p = ptp allows the com- 
mutation of the explicit factor p coming from the dc through the Taylor operator i2 + , that 
is assigned by eq. (||7|) to the 1PI part pln(p,m) of the above integral. For a = the new 
1PI integral is seen to be minimally subtracted -indeed, In(p,m) is logarithmically divergent 
(cf. appendix [A])- while for a = 2 it is oversubtracted. We are interested in the difference, see 
eq.( [5.25 ), i.e. in the subtraction terms themselves, of zero and second order in p: 



A2(z - y) ~ k {x - y) 



h 2 D 



d2 P Jv(x-y) P 2 



(2nf 



pp + m 



2 (tp - t 2 p ) ■ I u (p,m). 



(5.26) 



All terms are UV convergent, since they correspond to oversubtractions. The integral can be 
evaluated without further ado 



(t t) p -t 2 p )-I 11 (p,m) = -m 2 p 2 



d 2 



qq 



(2tt) 2 {qq + m 2 ) A 24vr m 2 
providing the following result for the difference ( |5.26j ) 



P 



A 2 (x - y) - A (x - y) 



ih 2 D f d 2 p 



r> 



24-7rm 2 J (2n) 2 pp + m 2 

After comparison with the rhs of ( |5.25 ) in the momentum representation, 



p{d 3 c(x)b(y)) cfi = -iph 



d2 P c ip-(x-y) P A 
(2tt) 2 pp + m 2 



(5.27) 
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the corresponding value for the contribution of the matter sector to the coefficient p in eq. ( [5.25 ) 
is read off: 

-KD 
Pm ~ 24W' 

In the same way, the ghost sector contribution may now be evaluated by considering the 
double contractions between T g ^(x) and @ g h(y) ( |5.23 ) in eq. (|5,25|) , taking the factor 
internal. The relevant BPHZ renormalized correlation functions A a (x — y) is 



K(x - y)=jr I d 2 z ( N 2 [2bdc + (db)c) (x)N 



b){(dc)c}+ (-c) {2bdc + (db)c} 



c,0 



-2h 7 



(27T) 



pp + m 



with the unrenormalized expression of the 1PI integral given by 



I(p, m) 



d 2 q 



2p 2 q + q 2 p - q A 



(2tt) 2 (qq + m 2 )[(p — q)(p — q) + m 2 } 
The ghost sector contribution to the lhs of eq. ( 5.25[ ) is the difference 

d2 P Jp-(x-y) P 2 r.l 



J^2{x - y) - AoO - y) 



-2V 



(2ny 



pp + m z 



(5.28) 



(5.29) 



Keeping in the Taylor series for the integrand in (5.28) only those numerators that survive the 
g-integration, namely the ones of the form (qq) n , yields 



(^ - tj) • I(p, 



m) 



d 2 q 



2qq 



(2tt)2 

The final result for ( |5.29| ) is then 

A 2 (x - y) - k (x - y) 



+ 



(QQ? 



(qq? 



(qq + m 2 ) 3 (qq + m 2 ) 4 (gg + m 2 ) 5 



13f p 3 
247r m 2 ' 



-26ih 2 



247rm 2 7 (2 pi) 2 



d2 P -ip-fa-H) 



pp + m 2 



from which, by comparison with the momentum representation ( |5.27 ) of the rhs of eq. ( |5.25D , 
the following value for the ghost sector contribution to the coefficient p in eq. ( |5.25| ) is found 



2QTi 



247rm 2 ' 

The results can be summarized by giving the Zimmerman identity ( p. 22 ) explicitly as 



N 2 [0(z)] - N [Q(x)\ = {2 \ A D)n — 2 ^2 



b*d A c(x) 



Inserted in eq. ( |5.21| ), it provides, after canceling the m 2 factors, the following form for the 
anomaly ( 5.20| ) 



A 



24vr 



d 2 x cd A b*, 



in full agreement with the result computed in the previous subsection (5.19). 

In conclusion, the local form of the anomaly equation enabled us to trace back the source 
of the anomaly to the oversubtraction of the variation of a mass term. The remainder of this 
evanescent term in the massless limit was extracted through a Zimmerman identity. 
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6 The anomaly of chiral W3 gravity 



In this section we first compute the anomaly of the chiral W3 gravity model [37] using the 
anomalous Ward identities ( [4.7D . Afterwards we also show how the two- loop contribution arises 
from (part of) the BRST breaking produced by the infrared regulating mass term, again through 
oversubtractions and Zimmerman's identity. 

VF-gravity models have recently attracted some attention for several reasons: as models 
with a non-linear symmetry algebra, as possible higher spin extensions of gravity, etc. For our 
purposes, the basic interest in the W3 model lies in the fact that it is an example of a field theory 
with a genuine new contribution to the anomaly at the two-loop level. In other words, the rhs of 
( |2.8[) contains at the two-loop level not only nonlocal one-loop dressings of the one loop anomaly, 
but extra local two-loop contributions to the anomaly. The expression of the two-loop anomaly 
on the local level, eq.( [4.9| ), rather than through its insertion in the effective action, adds extra 
clarity in that it is not necessary to disentangle it from the dressing of the one-loop anomaly. 
The expression we obtain for the anomaly in this model completely agrees with the existing 
literature [^, |39|, |2(], 41, 42, ^3], [5|, Existing derivations of the two-loop contributions 
always involved -implicitly or explicitly- the use of OPE's of conformal field theory (for the 
one- loop anomaly, derivations using traditional field theory methods are found in [pT| , ||). Our 
computation makes no difference between one and higher loop anomalies. 



6.1 The W 3 model 

Chiral W3 gravity is the minimal higher spin extension of chiral W% gravity, in which matter 
fields couple to gauge fields through their spin 2 energy-momentum tensor T and a spin 3 
current W 



1 



w 



1 



3 13 



d ijk (d^)(d^)(d<p k 



(6.1) 



The constants dijk that determine the W current are totally symmetric and satisfy 



for some arbitrary, but fixed parameter k. 

The extra gauge symmetry generated by the spin 3 current W is handled by supplementing 
the set of fields and antifields of chiral W2 gravity with an extra ghost-antighost pair (u, v), and 
their associated antifields. A tractable form of the extended actional for chiral W3 gravity in a 
gauge fixed basis can be taken from 00 



S 



d 2 x 



- (9^)09^) + b{dc) + v(du) 

+(f>* \c{d^) + ud ijk (d(f> j )(d(l) k ) - 2kb(du)u(dft 
+b* [-T + 2b{dc) + (db)c + 3v(du) + 2(dv)u] 
+v* [-W + 2kTb(du) + 2kd(Tbu) + 3v(dc) + (dv)c] 
+c* [{dc)c + 2kT{du)u] + u* [2{dc)u - c{du)}} 

S + J d 2 x <S>* A (x)R A (x) . 



(6.2) 



The antifields b* and v* take over the role of (minus) the gravitational gauge field and its spin 3 
counterpart, respectively. The expressions multiplying these antifields in (|6.2| ) are then (minus) 



20 See footnote |l| 

21 The free parameter a that appears in |2l| is taken to be zero here. 
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<S>* A R A 



the total energy momentum tensor Tt and the total spin 3 current Wt ■ The extended action 
shares an important property with (5.1), namely they are both of the form S = Sq+ J N2 
with So purely quadratic in the fields: correlation functions for this model can therefore again 
be considered as perturbative series in the antifields. 

The assignment of the relevant additive charges as in the W% model (Table ||), are now 
completed with Table for the extra fields. 





dim d 


spin j 


gh. (t 


d-j 


u 


-2 


-2 


1 





V 


3 


3 


-1 





* 

u 


4 


2 


-2 


2 


* 

V 


-1 


-3 





2 



Table 2: Additive charges for fields and antifields of the spin 3 sector. 



We finish the specification of the model by choosing the IR regulating mass term to be used 
in perturbative computations. For the matter fields the mass term (15.51) is used, such that the 



propagator is given by (5J3), whereas both for the spin 2 (b, c) ghosts and for the new spin 3 
(u, v) ghosts a mass term of the form ( fT7| ) is taken, leading to a propagator of the form (5.8) 
for both ghost pairs. 

6.2 W3 anomalous Ward identities 
6.2.1 Preliminary considerations 

We first analyze the anomalous Ward identities ( [4.7| ), as in the W2 case, to find the general 
structure of the W3 anomaly and describe the strategy we follow for its computation. 

As for the W2 case, it follows from Table [l] and Table [2] and the specific antifield dependence 
of the action that the BRST anomaly is linear in antifields, of the form A = <&\F A . The equation 
from which one can derive information on the function F A is again of the for m (|5.9P , in which 
the quantities R A should now be read off from the extended gauge-fixed action (|6.2D . Since loop 
diagrams come from at least double contractions between R (x) and R B (y) on the lhs of ( |5.9| ), 
the fact that R A contains terms with up to four fields indicates that relevant diagrams for the 
computation of the anomaly may contain up to three loops. However, a closer look at the 0(<I> 4 ) 
terms in the R A, s, given by 

-2kb{du)u{d(f> i ), 2k[Tbdu + d(Tbu)], 2k[T(du)u], (6.3) 

and at the form of the propagators, shows that no proper three-loop diagram can arise formed 
by quartic contractions between the expressions in fl6.3|) . The only two loop diagram arises from 



a triple contraction between two matter spin 3 currents W (6.1), so that the potential two- loop 
anomaly depends only on the antifield v* . Finally, double contractions between the different 
terms in the R A, s give rise to various one- loop contributions to the anomaly. Those are described 
in detail below. In summary, only one and two loop anomalies are possible in this model. 

The antifield coefficients F A (&) determining these one and two loop contributions to the 
anomaly are determined as in the W2 case. First one selects those pairs R A {x), R B {y) that 
admit double and triple contractions, and in a second step determines the "test product" of 
fields to be contracted with the remaining external legs. The precise form of the test prod- 
uct and a subsequent dimensional analysis completely determines the fields and the number of 
(^-derivatives present in the coefficient F A . The precise derivative structure and the numerical 
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coefficients are obtained from an explicit (renormalized) computation of the lhs of (|5.S|). In doing 
this last step, the minimal normal ordering of the interactions prescribes minimal subtractions 
for all the divergent one loop integrals, which should therefore be substituted by their mini- 
mally renomalized counterpart (|A.9|) . BPHZ renormalization of the relevant two- loop integral, 
requiring the use of Zimmerman's forest formula, will instead be treated in situ. 



6.2.2 One loop anomaly: a sample computation 



In this subsection, the procedure is illustrated by computing the "gravitational" part to the 
complete one loop BRST anomaly, viz. the contribution proportional to the antifield b* . This 
suffices to get a flavour of the method. Further details on the computation of the full one loop 
anomaly can be found in appendix |B|. 

Three different contributions to this gravitational part can be discerned: i) double contrac- 
tions of R b = —Tt with itself; ii) contributions containing one factor dij k ; and iii) contributions 
proportional to k, containing two such factors. We discuss them in detail in this order. 

i) This is analogous to the computation in section |5[ One uses R b = — (T + T2 + T3), with 
Tj defined by fl5,3| ) with suitable substitution of (b, c) by (u, v) in T3. Going through exactly the 
same steps as in sect. 5.2 , but taking two contributions of the type (|5.16| ) into account, namely 
for j = 2 and for j = 3, this first contribution to the one-loop anomaly of chiral W3 gravity is 



i(D - 100) 



24vr 



d 2 x cd 3 b* 



(6.4) 



ii) Contributions to the 6*-anomaly proportional to the tensor are generated by loops 
of matter fields coming from double contractions of — T in R b with the terms 

j ud ijk {d&){d(/) k ) in B* 
\ -id ijfc (^)(5^)W fc ) in R». 



Obviously, a suitable test product for both contributions is (j) l (y)v(z), so that this part of the 
anomaly has to contain one factor (j) 1 and one factor u, while dimensional arguments require four 
derivatives d. For this contribution, the Ward identity Q5.9|) takes the form 



- ( N 2 [-T(x)] N [ud ijk (d^)(d<p k )(y)\ v(z) 



-T(x)] N 3 

-«2 ( N2 



-djuid^id^id^z) 
F(^u,d A ;x)} f(y)v(z) 



c.Q 

(y) 



c.Q 



c,0 



(6.5) 



which can be represented by figure ||. These one-loop contributions are essentially the same 
as for the W%. They are minimally subtracted, i.e. up to quadratic terms in the external 
momentum (see eq.(|5.4|)). The resulting finite integrals are computed in eq.( [A.S| ), and the lhs 
of (16.51) becomes 
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Figure 2: One-loop anomalous diagrams for W3. The insertion in the first diagram is 0(y) = 
[ud ljk (d^){d^)]{y) 

This momentum dependence leads to the conclusion that in eq.( |6.5| ) F = d 3 (u(d(j)' l )diii) 1 and 
determines the coefficient 0*2- This yields the second contribution to the anomaly: 

Af ] = ^Jd 2 x [2u{dcl> k )d ijk ]dZ[b*8V] . (6.6) 

iii) Contributions to the 6*-anomaly proportional to k are generated by matter loops 
coming from two contractions of — T in R b with the terms 

2kTb(du) + 2kd(Tbu) in R v 
2kT(du)u in R c . 

A suitable test product for both contributions is now v(y)v(z)c(t), yielding a possible anomaly 
containing two u fields, one b field and four derivatives d in order to match the dimensions. ( |5.9| ) 
reads now 

{| ( N 2 [T(x)} N 3 [2kTb(du)(y) + 2kd(Tbu)(y)} v(z)c(t) ) cfi - (y ~ z)} 

+i ( N 2 [T(x)\ N_! [2kT(du)u(t)] v(y)v(z) ) c>0 (6.7) 

= -a 3 < N 2 [F(b, u, u, a 4 ; x)] v(y)v(z)c(t) ) cfi . 

The subtraction degrees computed from the -/V a , once more the minimal ones, namely 2, yield a 
renormalized expression for the lhs of ( |6.7l ): 

K 3 ikD f d 2 l f d 2 q f d 2 r i(q+l+r) . x j q . y ir . t u . z (g + I + r) 3 (g - I) 
12tt 7 (2vr) 2 y (27r) 2 y {2^ qlf 

1 
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The momentum dependence determines F = d 3 (budu). After computing the numerical coeffi- 
cient 03, the third contribution to the anomaly is 

Af i] = ^J^ X [-2fc6(9u)^']a 3 [6*5 ij ] . (6.8) 



In summary, collecting (p,4|), ( p,6|) and (|6.8|) , the gravitational part of the one-loop W3 
anomaly is 

Af = Af >m + < gh = ^ J d 2 x c^ 3 (6*«%) - ^ / d 2 ^ cd 3 6*, (6.9) 
where we have defined 

c ij = [c - 2kb(du)u]S ij + 2ud\{d(t) k ) (6.10) 

for further use. The subscripts m and gh will stand from now on for purely matter induced 
contributions and contributions induced by loops involving (one or two) ghost lines, respectively. 

This procedure is applied time and again in appendix |B| in order to find the complete one-loop 
BRST anomaly. The final result is 

At = Ai, m + < gh + Af* gh + Af jgh , (6.11) 



with w4^ gh , -4.1 g h and A[ gh given by (p\S|), (B.3) and (B/7) respectively, and where Ai, m lumps 



together all the matter induced contributions in the compact expression 

A 1 , m =^Jd 2 xcVd 3 h* j , 
with the "effective" antifield h*j given by 

ft* = {fy [b* + 2kb(u(dv*) - v*(du)) + 2kc*u(du)} - + 2v*d ijk (d(f> k )} . 

Our result for the one-loop anomaly in the W3 model is in perfect agreement with previous 
computations in the literature Q ^ |4l|, |2|, Q ||. 

6.2.3 Two-loop anomalous Ward identity 

After this brief discussion of the one loop anomaly, we now turn to the two loop anomaly. 
We consider it to be a major challenge for the present setup, since, as far as we know, no 
computation has been given that uses only the methods of renormalised perturbation theory. 

As already pointed out, the unique two-loop diagram constructed from three contractions 
between two terms from the available R A is the diagram that corresponds to a triple contraction 
between two copies of the spin-3 current W. This immediately leads to a two loop anomaly 
proportional to v* , since in the general Ward identity ( p\9| ) the i?'s are simply reduced to W. A 
suitable test product is hence v(y), which on the rhs yields a function F of the form ad 5 u, the 
only expression with non-vanishing correlation function with the test product v and with the 
correct dimension 2. The Ward identity ( |5.9[ ) takes the form 



* ( N 3 [W(x)\ N 3 [W(y)\ ) = a ( N 3 [d 5 u(x)] v(y) ) . (6.12) 



n 



This is the analogue of ( |5.10| ) in chiral Wi gravity and corresponds to the triple contraction in 
the standard conformal field theory OPE of W(x)W(y). 
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Figure 3: The anomalous two-loop diagram. 



The computation of the lhs of ( |6.12[ ), taking into account the appropriate number of subtrac- 
tions, is a straightforward application of the forest formula. It is given explicitly in eq.( |6.17| ). 
For readers unfamiliar with the rationale behind this forest, we build up the expression by con- 
sidering the regions in momentum space which cause a divergence in the unsubtracted diagram. 

The unregularized version of the two- loop integral E R (p,m) is given by 

I(p,m)= [ / ^-^t 7Z(ki,k 2 ;p,m), (6.13) 



(2vr) 2 7 (2vr 
with the integrand 

H(ki k 2 -p m) = - ^ - ^ (p - fci -k 2 ) 2 ^ 

k\k± + m 2 k 2 k 2 + m 2 [(p — k\ — k 2 ){p — k\ — k 2 ) + m 2 ] 

The momentum labeling is as in figure || In terms of the corresponding properly subtracted 
expression, the lhs of eq.(6.12) is 



~\ d2%2 i " -0f eiP ' {X ~ y)l[R (P> m )> with d 2 = d^ k d ijk . (6.15) 

Thus we have to determine that expression. We start by identifying the possible superficial 
divergences. There are four regions in the integration domain where the integral potentially 
diverges: 

k 2 finite quadratic divergence 

ki finite quadratic divergence , . 

k 2 — > oo, k\ + k 2 finite quadratic divergence 
k 2 — > oo quartic divergence. 

The quadratically divergent regions are associated with the one-loop subdiagrams ji, i = 1,2,3 
-pictorially represented in figure obtained by cutting one of the three internal lines between 
W^x) and W{y) in the original two-loop diagram (we take 71 to correspond to cutting the k 2 line, 
72 the k\ line, and 73 the (p—ki — k 2 ) line). The number of subtractions is, as always, determined 
from ( |5.4|) , and amounts to the minimal subtraction of the quadratic sub divergences identified 
above. The cut line corresponds to a field d(p, and if we include the resulting momentum 
dependence of the numerator (but not the propagator) in the expresssion for the subdiagram, the 
subtraction degree is 4. The quartic divergence of the complete diagram G also has subtraction 
degree 4. Zimmerman's forest formula [34] then consists of two steps. In the first step, the 




superficial quadratic divergences in the original integrand associated with the subdiagrams 7^ 
are subtracted: 



3 

4 = 1 
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Figure 4: The 3 one-loop divergent subgraphs. 



where p 7i stands for the external momenta of the ji subdiagram and 2T(7«) for the suitable 
restriction of the original integrand ( |6.14 ) to the subdiagram 74. This leads to an expression 
that is still quartically divergent. A further quartic subtraction cures this: 



(i - 4) 



1=1 

The net result of this process is Zimmerman's forest formula adapted to our integrand ( |6.14 ) 

(6.17) 



i=l i=l 

where T stands for the set of all possible forests IA of the two-loop diagram, namely 
T = {0, {G}, {71}, {72}, {73}, {71, G}, { 72 , G}, {73, G}} . 

Expression (6.17) is the subtracted integrand which upon substitution in ( |6.13 ) yields the proper 
expression of our two-loop integral. 

A closer inspection of ( |6.17 ) reveals that in its last two sums the terms corresponding to 
i = l,2 cancel each other. We exemplify this for 71, i.e. (1 — tp)tp H. The expression of the 
integrand ( |6.14|) restricted to the subdiagram 71 reads 



2T( 71 



(p + 77 



/c2^2 + fn 2 k\k\ + m 2 [{p + 77 — k%){p + 77 — k\) + m 2 ] 



where in the set of external momenta, p 7l = {p, ri}, 77 must be identified with — ^2 after having 
done the Taylor expansion. In this way, the contribution generated by the forest {71} is 



/1 9 



(6.18) 



where 0(p 2 ) stands for a polynomial of degree two in the momentum p whose coefficients are 
rational functions of k\ and k 2 . It is then evident that the action of (1 — tp) on ( |6.18| ) produces 
a vanishing result, since on quadratic polynomials in p, tp acts like the identity. The same 
mechanism can be seen to hold as well for the subdiagram 72, obtained from 71 upon interchange 
of k\ and k<i- 

For 73 the situation is slightly different. Indeed, the Taylor series expansion in the external 
momenta p 73 = {p, 7-3} of the integrand (6.14) adapted to 73, i.e. 



2T(73 



{p + r 3 - k 2 ) 



k2k2 + m 2 [(p + 7-3 — feXp + 7-3 — fo) + rn 2 ] [(p — k% — — h — fe) + m 2 ] 
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yields a nonpolynomial dependence on the external momentum p, namely 



t i Pi lL{^)\ T . A={ - p+kl+k2 ) 



k 2 



(p-ki- k 2 f 



k 2 k 2 + m 2 [(p — ki — k 2 )(p — k\ — k 2 ) + m 2 ] 



P13 



(p Ti -k 2 f 



[(P73 - k 2){Pyi -k 2 )+m 2 



p 73 =(fci+fc 2 ) 



This produces a nonvanishing result when the operator (1 — tt) acts on it. Using the more 
convenient integration variables I = k\ + k 2 and k = k 2 , the renormalized expression of the 



two-loop integral (|6.13[ ) is thus found to be 

d 2 l 



H R {p,m) 



(2vr)2 
d 2 k 



i-t 4 
2 {} ~ ) 



{p-if 



V > [{p-l){p-l)+m 2 



{I - kf 



kk + m 2 [(l-k)(l-k) + 



(6.19) 



It is now a straightforward task to evaluate ( 6.1S| ). Indeed, since the k— "one-loop" integral is 
precisely of the minimally subtracted form (A.1,A.2) considered in appendix [A|, with a = b = 2, 
use of ( A. 8 , A.9| ) brings it to the form 



f TL 



(p-if 



11 J 11 + m 2 [(p-l)(p-l)] + 



(6.20) 



Apart from the extra factor F (?p), with F(0) = 1, this is of the same generic minimally 
subtracted form (A.l. A.2| ), now with a = 4, b = 2. Therefore, in the limit of interest m — > 0, 
the one-loop result ( A.9[ ) can again be used, yielding for (6.20) 



1 p b 



480vr 2 p ' 

and for the lhs ( |6.15| ) of the Ward identity ( 6.12 ) the final result 



d 2 h 



d2 P c ip(x-y)P 5 

720vr 2 J (2tt) 2 p 
Comparing the expression ( |6.2l| ) with the rhs of ( |S.12 ) 



(6.21) 



iah 



(2vr) 2 



P 



the value of the coefficient a is read off, and the complete form of the W3 two-loop anomaly is 

A 2 = ^Jd"xud\\ (6.22) 

in agreement with |^, [|2|, ||]. This concludes our analysis of the one- and two-loop anomalies 
in the W3 gravity model by means of the anomalous Ward identities. 

6.3 The W3 two-loop anomaly from the mass term 



We show in this section how the two-loop contribution ( |6.22| ) to the BRST anomaly for W3 
arises as a consequence of the explicit BRST symmetry breaking by the IR regulating mass 
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term, i.e. from the eq. ( |5.20| ). The main tool is again the reduction of the oversubtracted BRST 
variation of the mass term through the use of Zimmerman identities. 



The W3 IR regulating mass term is formed by adding (|5.5D, (5.7) and a copy of the latter 
for the uv ghost system. Its BRST variation is 



N 9 



<5 S n 



5S 



5$ A (x) I 5$* A (x) 



-m 2 N 2 



<f(x) {c(d(f> 1 ) + ud ijk (dft)(d(f> k ) - 2kb(du)u(d<j> i )} 



+m^N 2 



g») {<?) + (gc) + 0) {«»} + Ou) {-W - W gh} 



(6.23) 



When inserted in loop correlation functions, it gives rise to a multitude of one- loop diagrams 
upon double contractions with the interaction terms <&* A R A in fl6.2|). These diagrams are of 
course the analogues of the one- loop correlation functions studied in sect. 6.2.2 and in appendix 
^. We will not dwell on a rederivation of the one loop contribution to the W3 anomaly, but 
concentrate instead on the more interesting two-loop contribution in this approach. Just as 
before, two-loop contributions can only be generated by triple matter contractions. These occur 
now between the interaction term — u*N3 [W] and the oversubtracted (f> 3 type terms in Q6.23 ), 
namely 



m 2 N 2 



f {ud ijk (d<tP)(d<i> k )} 



and 



m 2 N 2 



(6.24) 



Both contractions leave a factor u free and the two-loop anomaly is therefore necessarily of 
the form A 2 = m~ 2 pv*d 5 u, in accordance with the analysis presented in sect. 3.2.3| . Anisotropy 
forces the second term in ( |6.24D to behave effectively as a minimally subtracted, isotropic normal 
product, since the factor d~ 1 u is always in an external line. Upon substitution in eq. (5.20), no 
two-loop contribution is obtained from this term. The two-loop anomaly can then be extracted 
from the the Zimmerman identity 



N, 



^ {ud ijk (d<P)(dcf> k )} (x)] -N ud ijk {d<jP){dct> k ){x)\ =pN 



v*d 5 u 



+ 



(6.25) 



More specifically, the relevant coefficient will arise from a two-loop calculation. We define^] 



^a(x - y) 



d 2 z (N 3 [W{x)\ N a |y {ud ijk (d(fj)(dcp k )} (z)] v(y) 
The coefficient p in (6.25) is then determined by 

S 2 (x - y) - S (x - y) = p { N 3 \d 5 u(x)} v(y) ) . 

\ L J / c,0 

The unsubtracted 1PI two- loop integral corresponding to ( |6.26| ), is 

d 2 p 



2ih 3 d 2 



Jp-(x-y)_ 



where the integrand is 



H(p, m) = 
d 2 k 2 k 2 



(2vr) 2 

d 2 fci 
(2tt)2 



JP 

pp + m' 



: E(p, m) 



(6.26) 



(6.27) 



(6.28) 



H(k 1 ,k 2 ;p,m) 



A- 9 



{p-ki- k 2 ) 



(2tt) 2 J (2tt) 2 k\k\ + m 2 k 2 k 2 + m 2 [{p — k\ — k 2 )(p — ki — k 2 ) + m 2 ] 



(6.29) 



J Note that for minimally subtracted normal products, anisotropic and isotropic versions coincide. 
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which differs from the previous two-loop integral Q6.14 ) in one factor (p — k\ — hi). This changes 
the analysis of the divergent integration domains slightly in comparison with ( |6,16[) . Using the 
same momentum labeling, one has 




oo, k,2 finite 

oo, k\ finite 

oo, &2 - ► OO; hi + k2 finite 

oo, k2 — ► oo 



linear divergence 
linear divergence 
quadratic divergence 
cubic divergence. 



(6.30) 



The same labeling for the three different subdiagrams as in the previous subsection is used. In 
particular, 73 corresponds to the quadratic divergence in fl6.3C| ). The subtraction for the full 
two-loop diagram is always isotropic, and the degree is 2 + (3 — 2) + (a — 2) + (3 — 1) = 3 + a. For 
the subdiagrams 71 and 72, the (j) 1 is always in the loop itself, and the subtraction degrees are 
the same as for the two-loop diagram. For the third one-loop subdiagram (jf is an external line, 
and the prescribed subtraction is always the minimal subtraction to order 3, for both values of 
a. A straigtforward application of Zimmerman's forest formula on the integral fl6.2§ , 6\29| ) yields 
for the integrand determining the difference X2 — Sq Q6.27|) the symbolic expression 



8=1,2 



8=1,2 



JL. 



(6.31) 



The contributions related to 71 and 72, namely the last two terms in ( 6.31] ), are identically zero 
by exactly the same argument that was given in the previous subsection (cf. eq. fl6.1S| )). The 
difference in eq. (|6.27| ) is 



S 2 (x - y) - S (x -y) = 2ih 3 d 2 
with the 1PI part given by 
E R (p,m) = 



d 2 p 



Jp-(x-y)_ 



P 



pp - 



M R (p, to) 



(6.32) 



d 2 l 
(2vr) 2 
d 2 k 
(2vr) 2 



(tp tp^j 



ip-l) 



l-t; 



[ip-m 

kk + m 2 [(I - k)(l- k) + m 2 



I) + m 2 ] 
(I - k) 2 



(6.33) 



The same integration variables as in ( 6.19| ) were used. The consequence of the anisotropic 
normal product is that the subtraction of the one loop sub divergences is the minimal one. For 
the overall subtraction, there is a difference and the entire two-loop anomaly is due to this. 

We now compute Er. Using eq.( |A.8|) for the one-loop, minimally subtracted k subintegral, 
we obtain 



E R (p,m) = - — 



dx x 2 (l — xY 



d 2 l 



tt - 1: 



(p-i) 



[{p-l){p-l) + m 2 } 



x(l — x)ll + m 2 



Only the terms of the form {ll) n in the difference between the Taylor series survive the I angular 
integral, such that (|6.33| ) becomes 



Xr(p,™,) = —p 5 m 2 
Zn 



dx x 2 (l 



d 2 l 
(2vr) 2 



1 



x(l — x)ll + m 2 



(iiy 



(11 + m 2 ) 6 



The Wick rotated I momentum integral is conveniently performed using polar coordinates, and 
yields after changing radial integration variable to z = \l\ 2 /2m 2 , 



K R (p,m) 



1 



P 



4vr 2 



dxx (1 



dz 



(z+ l) 6 [x(l - x)z + l) 
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A further change of variables z = t/(l — t), followed by the change to 
{«i = tx, ct2 = t(l — x), 03 = 1 — t} brings it in the more familiar form 



F R (p,m) 



1 p b 
4tt 2 m 2 

1 p 5 1 

47T 2 m 2 



II / daAsiY.oci-l 



\i=l 



^(n/><)*(|>-i) 



aia2 + ai«3 + 0203 



where the last line is obtained by symmetrising the previous numerator in the on. Further use 
of the well-known result 



gives 



1 p 5 1 



27r 2 m 2 6! 

When plugged back in ( |6.27 , |6.32| ) this determines the coefficient p in the Zimmerman identity 

<i 2 ft 2 



&nm 2 



(6.34) 



Insertion of this Zimmerman identity in eq. ( 6.23|) leads then to the cancellation of its overall 
m 2 factor with the m~ 2 factor in (6.34) and finally gives a nonvanishing two-loop contribution 
to the anomaly ( 5.20| ) in the limit m — > 0. We find 

ihd 2 



. ihdr f l2 .t: . 

A * = m^J dxudv > 



which is in complete agreement with the previous result ( |6T2l ). 



7 Conclusions and Outlook 

In this paper, we have started a program to put the quantum theory of the BV quantisation 
formalism on a sounder basis. The old machinery of BPHZ renormalisation and normal ordered 
products turns out to be well suited for this purpose. We have given a formulation of the 
quantum BV master equation using these normal ordered products, which is valid to all orders 
in perturbation theory. An advantage over the Zinn-Justin formulation in terms of the effective 
action is that we work with local operator equations. As a concrete testing ground, we have 
computed the BRST anomalies for some two dimensional field theories, where it is known that, 
in contrast with Yang-Mills theories, new anomalies appear in higher loops. 

Apart from the possibility to compute the anomaly - to all orders in perturbation theory - 
directly from the Veltman-Ward identity, a closed local expression for the anomaly was obtained 
(eq.( [4.9[) ). Two ingredients are of crucial importance here. First, to apply the BPHZ method 
to massless theories one should temporarily introduce a mass term: the anomaly is related 
to the explicit breaking of the BRST symmetry by this mass term. This is reminiscent of, but 
different from, the one-loop Pauli-Villars method. The second ingredient is the use of anisotropic 
oversubtracted normal ordered operators. They originate in the mass term, and occur with a 
coefficient proportional to the mass (squared). The massless limit does not vanish however: it 
is extracted by rewriting the oversubtracted operator in terms of a set of minimally subtracted 
ones, using Zimmerman's normal product identity. The anomaly is what remains in this limit. 
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We have successfully applied our method to the computation of the anomalies in W2 and 
W3 gravity. In particular, detecting the presence of the two-loop anomaly of the latter model 
is as easy as detecting the one-loop anomaly of either model. Only the computation of the 
coefficient, involving the renormalisation of a two-loop diagram using forest formula techniques, 
is more involved. 

Of course, these methods can also be applied to compute anomalies in global symmetries. As 
an example, we verified that with our choice of mass term, the correct ghost number anomaly 
is reproduced. Despite the fact that both the original action and the mass term have ghost 
number zero, the associated Noether current is not conserved, owing to a non-locality in the 
mass term, and again the occurrence of oversubtracted operators. 

Having demonstrated the possibilities offered by the BPHZ method in setting up proper 
renormalised versions of the quantum master equation, we now give some possible directions for 
further development. 

A straightforward extension would be to include extra finite "counterterms" in the action. In 
specific examples, these harbour the possibility to cancel or shift the anomalies of the quantum 
theory. These are related to the possibility that the interaction terms and/or the BRST trans- 
formation rules are not minimally subtracted, and presumably also to the choice of mass term, 
as in Pauli-Villars regularisation. Also, whereas our local anomaly equation is in principle valid 
irrespective of the closing of the algebra off-shell, it could be of interest to work out explicitly 
an example where the extended action is non-linear in the antifields. 

On a more theoretical level the question of consistency conditions for higher loop anomalies 
imposes itself. For the effective action we have that (T,A ■ T) = . Formally, it is easy to derive 
the corresponding local conditions perturbatively in h, but unfortunately this is again ill-defined 
due to the use of the formal operator A (see section ||). To lowest order, this corresponds to the 
condition of Wess and Zumino that the one-loop anomalies are consistent, i.e. BRST closed. In 
the BPHZ formulation that we gave, the offending A is absent^. The question of higher-loop 
Wess -Zumino condition is therefore open. Our local anomaly formula, eq,(f4.S|), suggests that 
the normal product equations could be rather simple transcriptions of the classical equations, 
provided one uses the appropriate anisotropic versions. 

The definition of the normal products takes as its starting point a set of well-defined Feyn- 
man rules. This means, in BV, that one starts from a certain gauge fixed basis. It would be 
interesting to find a formulation that does not depend on this choice, or at least to show how 
the transition from one choice to another one can be made. This poses of course the problem 
of treating canonical transformations (in the antibracket sense) on the quantum level with the 
BPHZ method. We recall that within the BV scheme, a change of gauge is implemented via a 
canonical transformation. Possible counterterms that are generated by this are also related to 
the Jacobian of the BRST transformation. Although for a closed algebra it is formally indepen- 
dent of the gauge, in a regularised computational one may obtain a gauge dependent expression, 
owing to the gauge dependence of the regulator. 

The final task would be the unraveling of the quantum cohomology. On the non-local level, 
the quantum BRST transformation can be defined as (r, X) = SX. For non-anomalous theories, 
this is a nilpotent operator. It can be translated back formally to the level of local field theory, 
where the quantum BRST operator is aX = (X, W) — ihAX, again involving the ill-defined A. 
One may hope that the normal products of Zimmerman will enable one to define a replacement 
for this operation, not just for the action (this is what the anomaly gives), but also for arbitrary 
functionals of fields and antifields. 

23 Another possibility is to stay with A, but de- localise all the functionals, as in the nonlocal regularisation 
scheme treated in ||. See Q for further developments in this direction. 

24 We refer to the case of Yang-Mills regularised with Pauli-Villars, see |4J]. 
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A BPHZ renormalized one-loop integrals 

Our main goal in this appendix is to compute the general form of the BPHZ renormalized, 
minimally subtracted, one-loop integrals used throughout the paper. 

At the one-loop level, the specific form of the interactions for chiral W2 and W3 gravity leads 
to connected Green's functions which always contain as a 1PI part a special case of the general 
one-loop integral I a b(p,m) ( |5.13j ) 

r ( , f d 2 q q a {p - q) b _ f d 2 q 

Uip, m)=J-^y 2 {q _ - m2)[{p _ q){p _ ¥) +m2] =J J^y^P, m), (A.l) 

with a, b > 1. With a superficial degree of divergence given by a + b — 2, minimal subtraction 
prescribes that these integrals are replaced by 



I*(p, m) = [(1 - t a + b - 2 ) ■ I ab (p, m)] , (A.2) 



using the shorthand ( |5.12| ). 



First one expresses the propagators in ( A, 2 ) by means of their well-known integral represen- 

'CX 





tation a 1 = J °° alte at , i.e. 



After changing integration variables to t = t± + ti and x = tilt, this becomes 

l3,(p,m) = J ^2 / o °°d^y 1 dx e -^ +m2 )(l-^ +6 - 2 ) [q a (p-q) b e- tx ^-"^] . (A.3) 

For the actual computation of this integral it is not necessary to know the subtraction terms 
in detail, the knowledge of their general form is sufficient. The Taylor series expansion around 
vanishing external momenta of the relevant factor of the integrand in ( A.3| ) is 

q a (j) — q} b e - tx [pp-gp-gp] = 

b 00 r+s<n / h\ 1 

(_-t\k+r+s+nl 1 + n rr n „a+k+s =n-r b—k+n—s 

hhrh [k)rls\(n-r-s)\ tx ^ P P 

The terms that must be subtracted in the minimal scheme are determined by the inequality 

b-k + 2n-r-s<a + b-2, 

from which, upon use of r + s < n, the inequality a + k + s > r follows. All the subtraction 
terms have therefore the following q dependence 

q y (qqf, with y > 0, (A.4) 
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i.e. they always contain extra factors q, not paired with q. The renormalized, minimally sub- 
tracted integral (|A.3|) becomes then 

/«(p,m)= / 7^ [°° dtt f 1 dx e -*(^+™ 2 ) La( p _ 9 )6 e -tebP-#-w] -Y'a rs q y (qq) x p r p 
J [2ir) z Jo Jo 1 ' 

(A.5) 

where 2~2' stands for the sum restricted to the terms satisfying the conditions r + s<a + b — 2, 
y > 0. Treated separately, the terms in the integrand of ( |A.5| ) may be UV divergent, but 
keeping them together the integral is by construction UV convergent. The mass m ensures 
IR convergence (the region t — > oo) as well. In conclusion, ( |A.5| ) is an absolutely convergent 
multiple integral on which manipulations like changes of variables, interchanges in the order 
of the integration, etc, are mathematically well-defined, provided the integrand is treated as a 
whole. 

Now we use this freedom, and first perform the integral over q. Due to the general form 
( [A.4D of the subtraction terms and the general result 

|^-«V = 0, (A.6) 

none of them survives the q integration, yielding 

I*(p, m ) = £ dx J™ dtt J ^2Q a iP ~ q) b e -^ q - xp ^- x ^ +M ^ , (A.7) 

with M(x) = m 2 + x(l — x)pp and where now the order of integration is fixed. The net effect of 
the subtraction terms has been that they have imposed a fixed integration order. Changing the 
loop momentum integration variable to k = q — xp and using ( |A.6| ) , brings ( |A.7| ) in the simpler 
form 

rl roc f A^V, 

I*(p,m) = J o dx J Q dtt j ^e~ tkk e- tM ^x a {l-x) b p a+ \ 
which after trivial (Wick rotated) k and t integration (in this order) results in 

'Sfo m) = ^V+ b ^ dx 2 X " {1 ~ X) \ _ = I* 6 (P, 0) F ab (?C\ , (A.8) 
2tt Jo + x{l — x)pp \PP J 

with .F a &(0) = 1. The massless limit provides the desired result 

^S(P,0) = ^5_/ tear-^l-x)*- 1 ^ _L_B(a,6), (A.9) 
2tt p Jo 2-k p 

in terms of the well-known Euler Beta function B(a, b). Expression (A.9) is the general expression 
for the minimally subtracted one-loop integrals needed in our computations. 



B Appendix: The complete one— loop anomaly in W3 

In this appendix, the computation of the W3 one-loop anomaly is completed, by analyzing 
all the possible double contractions between R (x) and R B (y) in eq. (|5.9|) . The analysis is 
performed antifield by antifield, where the antifield in question determines the R A {x) to be used 
on the lhs of ( ^.9[) and the "$^-anomaly" or coefficient F A of the one-loop anomaly that is 



determined. The 6*-anomaly is computed in sec. 6.2. 
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u* anomaly 

The BRST transformation for u, R u = [2(dc)u — c(du)], can not be contracted twice with 
any other quantity in the R s: no terms with both a b and a v antighost are contained in them. 
Consequently, no term proportional to u* arises in the anomaly. 



c* anomaly 

The only term in the BRST transformation R c of c that can give double contractions with 
other terms in R B is the monomial 2kT(du)u. Moreover, this always goes via the matter loops 
studied for the 6*-anomaly in the main text. The direct consequence of this simple fact is that 
all the Ward identities that hold for the matter-loop induced part of the 6*-anomaly will hold 
as well upon replacing the antifield b* with the combination —2kc*(du)u. In fact, this amounts 
to rewriting the relevant interaction term c*[2kT(du)u] in ( |6.2| ) as the product of an "effective" 
antifield h* = —2kc*(du)u, with exactly the same quantum numbers as b*, times minus the 
matter energy momentum tensor T. It is then clear that the c*-anomaly is, in analogy with 

(PD 

A( jm = ^Jd 2 x c ij d 3 [2kc*u(du)5 ij ] , 
with defined in (6.10). The same argument allows a simple determination of the contributions 



to the 0* and v* anomalies induced by matter loops as well. 



<j)*— anomaly 

In the BRST transformation for the matter fields 

Bf = ud i jk (d(/> i )(d(l) k ) + c{df) - 2kb(du)u{d4> i ) , (B.l) 

the first monomial, can only produce contributions via matter loops, in analogy with the c*- 
anomaly. Substituting b*5{j in ( |6.9D by the effective antifield —2u<p k dijk leads then to the fol- 
lowing matter induced contribution to the 0* anomaly 

The remaining two terms in (|B.1| ) generate contributions to the tfi* anomaly through mixed 
matter-ghost loops or by purely ghost loops. Indeed, the monomial c(d(ft l ) gives rise to mixed 
matter-ghost loops through contractions with the terms 

-2/c6(9u)u(a0 i ) in R* 
2kTb(du) + 2kd{Tbu) in R v , 

while the monomial —2kb(du)u(d(p t ) produces also mixed matter-ghost loops through contrac- 
tions with the term c(dcf) 1 ) in R& and purely ghost loops through contractions with the term 
3v(dc) + (dv)c in R v . A good test product of fields for all these four possibilities is <fP (y)v(z)v(t), 
which forces the corresponding anomaly to contain one </> field, two u fields and four derivatives d. 
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The one-loop anomalous Ward identity Q5.9| ) takes the form 

| ( N [c(0^)(a;)] N [-2AA(0u)u(0^')(y)] v(z)v(t) > c0 

-| ( No [c(a^)(x)] N 3 [2fcT6(<9n)(z) + 2fc0(T&u)(z)] w(t)^'(j/) > c0 - (z, t) 

+i < N [-2M>(duMd^)(x)] N [c(^')(y)] u( Z )t;(t) ) cfi 

-| < N [-2tt(0n)n(^)(x)] N 3 [3«(9c)(z) + (0u)c(z)] ) Cj0 - (s, *) 

= -a 5 < N u, it, a 4 ; x)] & (y)v(z)v(t) ) c>0 . 

After a lengthy but straightforward calculation, the lhs is found to be 



7T 



d 2 l 



d 2 q 



(2vr) 2 J (2vr) 2 J (2vr) 



d 2 r 



e i(g+i+r)-x e -ii}-y e -ir-t e -ii-2 



-i(r - /)(/ + 9 + r) 2 - ~(r - /)c? 2 + 2(r - i)W + ^(/ 3 



(B.2) 



where the factor 1/g corresponds to a propagator between 50 and </> and the factors 1/f and 
1/1 are generated by (u,v) ghost propagators. The combinations of momenta in (B.2) can be 
associated to the following monomials 



I 3 — r 3 
r 2 l — l 2 r 
(r - l)q 2 
(r-l)(l + q + r) 



u(d 3 u)(d4>) 

(du)(d 2 u){dcj>) 

u(du)(d 3 (j)) 

d 2 (u(du)(d<j>)) . 



After computing the coefficients of the different terms in the usual way and performing a few 
integrations by parts, the second contribution to the eft* anomaly is 



a4> 



v* anomaly 



—ik 
6tt 



d 2 x 4>* 6d(u(<9u)(dy )) + V{d' z u)(du){d(t) 1 ) + &u{d 6 u)(d(t?) 



(B.3) 



Let us finally compute the contribution to the one-loop anomaly induced by the BRST 
transformation of the ghost field v 



R v = -W + 2kTb(du) + 2kd(Tbu) + 3v(dc) + (dv)c . 



(B.4) 



The contributions caused by matter loops are generated by double contractions between matter 
fields in —W or in the energy-momentum tensor T in (B.4) with (d(j)) 2 type terms in the R B, s. 
When writing the interaction terms in the form of effective sources times — T in analogy with the 
matter induced 6*-anomaly, the following matter induced contribution the one-loop u*-anomaly 
is inferred 



2v*d ijk (d(f) k ) + 2kb(udv* - v*du)5,, 



The computation of the contribution of the loops with at least one ghost line, requires more 
effort. Closed loops of one matter field and one ghost or of two ghosts can be obtained in the 
following four ways. The terms [2kTb(du) + 2kd(Tbu)] get contracted twice with 

f 3v(dc) + (dv)c in R v 
\ c{d(f) in R*, 
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while the terms [3v(dc) + (dv)c] can contract with 



j —2kb(du)u(d(j) t ) in 

| 2kTb(du) + 2kd(Tbu) in . 

A suitable testproduct for all four contributions is <ft (y)4> k (z)v(t), yielding a contribution to the 
v*-anomaly constructed out of two 0-flelds, one u ghost and five derivatives d, and a one-loop 
anomalous Ward identity of the form 



i ( N 3 [2kTb(du)(x) + 2kd(Tbu)(x)} N 3 [3v(dc)(t) + (dv)c(t)] $ (y)(j) k (z) 



c.Q 



+ f ( N 3 [2kTb(8u)(x) + 2kd(Tbu)(x)] N [#)(]/)] ^M*) ^ + (j, fc)(y, z) 



+| ( N 3 [3v(0c)(a:) + (&>)c(x)] N [-2kb(du)u(d^)(y)] ^ k (z)v(t) ) ^ q + (j, fe)(y, z) (B.5) 



jr ( N 3 [3u(5c)(x) + (dv)c(x)} N 3 [2fcT&(8tt) (i) + 2A;0(T6u)(t)] ft(y)(j> k (z) 



c.Q 



ae, 7 ( N 3 [F(u, <f>, 0, 9 5 ; x)] <j>>(y)4> k (z)v(t) 



c,0 



The equality of the first and fourth contribution on the lhs of ( |B.5| ) upon exchanging x and 
i, and the fact that the value of the second and third contributions can be read off from some 
of the results obtained in the (j)*- computation above, reduce the amount of work necessary to 
obtain the renormalized expresion for the lhs of ( B.5 ). It reads 

-kh 3 f d 2 l f d 2 q f d 2 r 



TT 



(2ir) 2 J (2tt) 2 J (2vr) 
1 

qlf 



e i{q+l+r)-x e -iq-y e -ir-t e -il-z 



Ur(l 2 + q 2 ) + (l + q + r)(l 2 + q 2 



+ - (5r 3 + \2r 2 (l + q + r) + 12r(Z + q + r) 2 + 5(1 + q + r) 
6 V 



(B.6) 



It is now a piece of cake to determine the distribution of the five derivatives d in the coefficients 
F(u, 4>, (j), 9 5 ; x). Indeed, taking into account that the factor (qlf) -1 is generated by two prop- 
agators between <f> and deft, and one propagator (u,v), the other momentum functions can be 
associated with the following terms in F(u, (j), <f>, d 5 ; x) 



-> (d(j))(d 3 (j))du 

-> d(u(d<p)(d 3 <t))) 

-+ (d 3 u)(d(j))(d<f) 

-+ d ((d 2 u)(d(i>)(d<l))) 

-> 5 2 ((an)(a0)(a</»)) 

-> a 3 (u{d<l>)(d4>)) . 

By simple comparison with (|B.6| ) , the coefficients of each of these terms are determined and the 
mixed and ghost induced contributions to the one-loop v*-anomaly is 



f r(l 2 + q 2 ) 

(r + I + q)(l 2 + q 2 ) 
j. 3 

r 2 (r + 1 + q) 
r(r + I + q) 2 
, (r + / + <?) 3 



■A? 



ik 



— / d 2 x (v*(0u) 

27T 



07T 



5(<9 3 u)t;* - 12(<9 2 u)(dv*) + 12(0u)(0V) - 5u(dV 



(B.7) 



39 



In summary, collecting all the previous contributions, we can write the complete one- loop 
anomaly for W3 in the compact way ( |6.11 ). This finishes our analysis of the one-loop contribution 
to the BRST anomaly in the W3 model. 



C The ghost number anomaly 



In this appendix, we show how the well-known ghost number anomaly [ 45 1 is obtained in the 
context of the BPHZ scheme. Although in principle this is another application of the method 
sketched in sect. |3.4| for computing the global axial anomaly, we find it instructive to present 
it as an example of how nonlocal mass terms, though invariant, may give rise to anomalous 
conservation laws. Also, the fact that the correct ghost number anomaly is obtained should 
dispell any misgivings about the use of this non-local mass-term. 

In order to make contact with standard formulations of chiral gravities, we consider a sim- 
plified version of the models used before, with all antifields set to zero except the antifield b*, 
which takes over the role of (minus) the gravitational gauge field, i.e. 



S 



d 2 x 



1 



m 
~2~ 



b* [T + T 2 + T 3 ] 



(C.l) 



{cj,bj}, j = 2, 3 are the spin 2 and 3 ghost-antighosts systems and T2, T3 their corresponding 



energy momentum tensors (|5.3|) . As it stands, this massive action (|C.1[) is invariant under the 
following ghost number transformations 



Sib* = 0, 



(C.2) 



where the ghost numbers $a for the fields are listed in tables [l] and One would naively conclude 
that there is no room for a ghost number anomaly in our formulation, contrary to the well-known 
results. However, this is not the case and, once again, the oversubtractions performed on the IR 
regulating mass term, together with its nonlocal character, are the origin of the ghost number 
anomaly. 

As usual for global symmetries, the analysis starts by computing the divergence of the ghost 
number current (S = S{m = 0) 



5S 



-5^ A (x) 



+ m 2 ^2 

.7=2,3 



\ b i) C J + b i {\ C 3 



5$ A (x) fl 

which classically vanishes on-shell in the limit m — > 0. The anomalous Ward identity is 



(C.3) 



N 2 [d^(x)}]l^(x t ) 



ffiEM n*^) ) s(x-xi)+( N 2 [g(x)]n^ 



(C.4) 



where the operator insertion Q(x) parametrises the ghost number anomaly. The linearity of 



the transformations ( C.2 ) allows, upon application of the Schwinger-Dyson equations (|3.5| ), to 
rewrite the contribution from the massless piece in (|C.3|) as 



N9 



5S 



S^ A (x) 
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which reproduces the first, non-anomalous term on the rhs of the Ward identity ( |C.4j ). The 
ghost number anomaly is therefore determined by the over subtracted operator 



m 



E N 2[(^)c,(x) + 6 i Qc,(x)) 



3=2,3 



No anisotropic normal orderings appear because the transformation rules of the ghost number 
symmetry are linear. 

Once again, prior to taking the limit m — > 0, the Zimmerman identity ( |3.7| ) should be used 



m 2 \ N 2 

3=2,3 



l b 3 ) O + 



1 

9 Cj 



N n 



Cj + bj 



1 

d Cj 



}=E r fe N 2[^]- ( C - 5 ) 

J I. 



Its rhs is determined in the usual way by inserting both sides of this relation in suitable cor- 
relation functions. In doing this computation, however, some work can be saved by using the 
dimensional arguments presented in sect. 5^. They show that the insertions Qk(x) in ( C.5 ) must 
be of the form Q/- ~ G(d)b* . In this way, it is concluded that the one-loop contributions to the 
correlation functions 



M>(x) = ~ I d 2 yb*(y)P a (x-y) 



- I d z yb*(y) ( N 2 \jbjidcj) + (j - l)(dbj) Cj ] (y)N„ 



-bj ) Cj + bj [ -Cj 



d 



c,0 



given in terms of the BPHZ renormalized one-loop integrals 



I3 a (x-y) 



d 2 p 



Jp-(x-y) 



(2vr)2 l L p' 



PQ + (j ~ l)p 2 



(qq + m 2 )[(p — q){p — q) + m 2 



(C.6) 



provide the relevant information on what the rhs of ( |C.5| ) exactly is in terms of the difference 
V- 2 JMU.r; M^x)). 

Computing the difference P 2 (x — y) — Iq (x — y) requires the knowledge of the p expansion of 
the integrand in ( |C.6|) up to second order. Keeping only terms of the form (qq) n in the Taylor 
series (see appendix |A]), we obtain 



d 2 p 
(2tt) 2 



Jp-(x-y) 



d2 » e 



{21,) 



2 V 



PQ + (j - 1> 2 



(qq + m 2 )[(p - q)(p - q) + m 2 ] 



d P Jp-(*-y) „2 
(2vr) 2 P 



d 2 q 
(2vr)2 



u - 1) 

(qq + m 2 ) 2 



qq 



(qq + m 2 ) 3 



K2j-m 2 f d 2 p M *-v) p z = jw-^ 2 dH{x _ y) . 



Airm 2 



(2n)< 



Airm 2 



Plugging this result back in the difference Sj=2,3(-^2 ( x ) ~~ -^o( x ))' we easily obtain 



M 3 2 (x)-Mi(x) 



(2j - l)h 
Anm 2 



d 2 b 



The Zimmerman identity (|C.5|) then takes the form 



m 2 £ \ N 2 

3=2,3 



~) h l I ' 



' > + h J [ \ C 3 



Nn 



Cj + bj 



(3 + 5)ft 2 
4tt ' 
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from which, upon taking the limit m — > 0, the expression of the anomalous insertion in ( |C.4| ) is 
finally determined 

g( x ) = —d 2 b*(x). 

TV 

This is in agreement with previous computations in the literature [fl5[ |. 

We conclude that, despite the invariance of the mass term, the correct anomalous Ward 
identity is obtained. The anomaly arises from oversubtracted operators which occur with a 
coefficient proportional to the mass (squared). The massless limit of these combinations is 
extracted using Zimmerman's normal product identity. 
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